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PREFACE TO FIRST EDITION 

The view is gaming ground in the educational world that to free the 
learner (after a certain point) from laborious calculation is to help him to 
grasp and develop the underlying ideas. These exercises are therefore 
devised to use the machine so as to stimulate mathematical thought. This 
aim has been pursued from the beginning but only becomes fully 
apparent after the full range of machine operations has been mastered. 

This collection of exercises, along with the necessary minimum of 
explanation, was devised for the use of sixth-formers. It should be 
helpful to ^osc teachers who want their pupils to be introduced to 
machines but who are not clear how to set about it. 

In the first five chapters the technique of handling a machine is 
developed stage by stage through progressive examples. For each stage 
it is necessary first to learn certain manual operations, preferably by 
seeing them demonstrated and immediately imitating them. For a 
course in being, one pupil can demonstrate to another — to their 
mutual advantage. A teacher wishing to initiate such a course can acquire 
the skills with aid of {a) a fiill demonstration by the maker’s representa¬ 
tive, followed by (i) thorough practice, with the maker’s handbook for 
reference. 

The author wishes to express gratitude, for their encouragement and 
help, to Mr, M. Shoenberg (educational adviser of Olympia Business 
Machines, who make the Brunsviga), to Mr, F. W. Russell (education 
manager, Monroe Calculating Machines), to Dr. J. Crank (Head of 
Mathematics Department, Brunei College), to Mr, M, N, Horsman of 
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PREFACB TO FIRST EDITION 


Brunei College, as well as to Mr, D, Grisewood of Macmillan & Co, Ltd, 
and to Mr, M, Bridger of Leicester College of Technology. 

Acknowledgement must also be made to the author’s pupils at St, 
Paurs School, who have greeted machines with enthusiasm, treated 
them with respect, and made shrewd comments on the exercises. 

Standing orders for pupils doing the course are reproduced as an 
appendix, together with some hints which may prove useful to teachers 
who wish to use a machine occasionally for junior class-teaching. 


A. J. M, 
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INTRODUCTION 


The Desk Calctilating Machine 

Machines operated entirely by hand invariably incorporate the 
following registers, where the numbers which are involved in the cal¬ 
culation are exhibited: 

(1) The input or s^tUng mghter, 

A number is set on the machine by the pressing of keys or the 
setting of levers; and the setting register indicates this number, as a 
check that (a) its digits have been set correctly and (5) it is correctly 
located* 

(2) The i^cle register or multiplier register* 

The processing of the input number is done entirely by turns of 
the operating handle, c*g. a single fonv^rd turn of this handle causes 
the input number N to be added in to the accumulator once (while 
still remaining set for further use)* This turn is recorded as a unit 
figure in the cycle register* Seven such turns would show the number 
7 and give 7N in the accumulator: this fact explains the name 
^multiplier register* often employed. 

(3) The acmmulator^ already mentioned, can also be called the pro¬ 
duct register as the above example explains; but the former term 
explains more clearly its usual function, viz*: to accumulate (so long 
as no action is taken to clear it) all numbers which, singly or as 
multiples, are brought into it by the action of the operating handle* 


The Operation of the Machine 

This is described in general terms which apply to any machine except 
where special reference is made to a particular model. 

The basic function of the operating handle has been described; but its 
operation would be very limited without the shift lever^ a device which 
by lateral movement of the accumulator causes the numbers to be fed 
into it at any position. E.g. with the input one place to the left (relatively) 
the operation described in (2) would feed in not 7Nhut 7QN^ 
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INTRODUCTION^—THE DESK MACHINE 


It must be possible to clear any register, and we have therefore three 
separate clearing leversj and sometimes a total clearance lever which clears 
all registers in one motion* (It is however possible, e*g. in the Bnmsvtga 
20, to lock the left-hand half of the accumulator as a storage, and the 
clearance will then not operate on this half)* 

Another valuable facility on many machines is a hack-transfer device by 
which a result can, while being cleared from the accumulator, also be 
automatically set. By adjustment of the shift, this setting can either be 
done complete or with the final digits omitted* Generally, all but the 
smallest machines now have this back-transfer facility. The possession 
of a machine with such a device eliminates one source of error (viz. 
resetting by hand) as well as speeding up certain programmes of 
calculation. 

A list of firms marketing hand-machines in the U.K* is given on 
p, 61, with some notes on choosing a machine for teaching purposes. 

Electrically driven machines are basically the same as the hand- 
machine in their mechanical construction, but fall into three classes 
according to their degree of automatic operation: 

(1) Semi-automatic machines, in which division is automatic but not 
multiplication; 

(2) fully automatic machines, in which both division and multiplica¬ 
tion can be preset j and 

(3) machines like the Marchant in which the multiplicand is set on 
the keyboard and the multiplication process is then carried out 
digit by digit. 

The possibility of jamming a machine is increased by automatic 
operation. It usually occurs because a new process is initiated before the 
machine has completed the previous one. 


INTRODUCTION 
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Note for the Learner 

To the learner; the course is planned as follows: 

Before each of the early chapters you must have certain machine 
operations shown to you on your own type of machine. These are listed 
in the notes at the appropriate stages. Generally, anyone who has suc¬ 
cessfully completed the exercises can act as demonstrator for you. 

Repeat the manual operations until you understand them and do not 
let the demonstrator go until you are satisfied. You should then do the 
exm'dse without help even if it means trying it a number of times. 

There are no ‘answers at the back’. In every case you will be able to 
verify your own numerical answers by the checking method indicated. 

Keep a special ‘maths practical book* for this work. Record all your 
results, both intermediate and final; together with your method if any 
choice was given, and your check. A clear and S 3 rstematic layout is of 
great importance, and in the more complex examples a specific form of 
layout is given for you to follow. 

Where there are verbal questions, answer them in the book. Add any 
useful comments of your own at the time: this will increase the value of 
the record. Solution notes are given on pp. 62 onwards, to questions 
which seem to demand them. The existence of such a note is indicated 
in each case thus; [Solution note. 

In professional numerical work tio answer is given without an indica¬ 
tion of its acceptability — by rounding off figures or otherwise. In this 
course you will be gradually subjecting your results to criticism; but you 
should from the beginning avoid copying down figures from the machine 
*Ju$t because they are there*. 

The book consists of 10 chapters. An o^ercise numbered 3.6 is the 6th 
exercise in the 8th chapter. The early explanatory notes precede the 
exercises to which they refer, but from chapter 4 onwards the left-hand 
pages are reserved for notes and the right-hand for the exercises. 

Considerable cross-reference is used where techniques arise again in a 
new context, but for additional guidance you will find on each page a 
short summary of the contents, and at the end of the book is a short 
index of terms, giving in each case the page where the term first appears. 
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CHAPTER 1 


SUMMATION AND THE DIFFERENCE TABLE 

Learn how 

(i) to dear all registers; 

(U) to return the carriage to first position ; 

(iii) to set decimal point markers for addition/subtraction; 

(iv) to set up a number on the levers or keyboard, and to verify that it 
is correctly set; 

(v) to add this number into the accumulator; 

(vi) to notice the indication given that an addition has been carried 
out; 

(vii) to dear the keyboard only, ready to set and add a further num¬ 
ber. 

LI* Add 592 and 89, in this order and note the result* Then add in the 
reverse order, 

1.2. Add 8-649 and 14*7 in this order. Note. 

Add 14'700 and 8*649 in this order. 

If the results do not agree, probably your decimal point setting was not 
correct : the decimal point marker, set for the largest number of decimal 
places (3 here), applies to all numbers set. 

1.3. Add each row of the following table, noting the results and then 
adding them to give a grand total. 

27*94 6*79 18-312 43*0 =sum of row 

31-02 28*4 43-89 174 « „ „ „ 

6-7S 5-8! 1*08 25-02= „ „ „ 

total 


Now repeat the summation, but add first the columns, and then sum the 
column totals to a grand total. 

What is the significance of the number which appears in the cycle 
register at any stage of the calculation 1.3 ? 








2 SUMMATION WITH SIGNED NUMBERS 

For 1.4, Learn 

(i) bow to subtract a number from the accumulator; 

(ii) what happens If the number subtracted is greater than that in the 
accumulator; 

(iii) when it Is possible to disregard the appearance of such results 
(see note on complemmis below). 

(iv) how to restate a number obtained as a complement, in a form 
with a negative sign. (Calculate mentally; and add back, as a 
check, to give zero.) 

1.4. Add the rows and columns of the table below in the same way as in 
L3, and note all results/Wi os th^ appear m the machine. Form the 
grand totals as before. 

72*94 --6*97 18-321 -40-3 

-13-02 28-4 43-89 -17-4 

6-75 -5-81 -1-08 25-02 

1.5. Convert the result of adding the 4th column of 1.4 into a negatively 
signed number. 

1.6. Set a number which includes a decimal point. Add it into the ac^ 
cumulator 11 times. 

(You have thus learnt how to multiply by an integer). 

Would it be possible to get this result by only two turns of the handle 
instead of 11 ? [Solution note. 

A note on complements 

If you subtract unity from an 8-digit (cleared) register it shows 
99,999,999; while 0-JV shows (10®-N) say. N' is called the 
complement of iV, since they make up together a standard total (cf. com¬ 
plementary angle in geometry). 

An all-machine method of finding N from N* (standing in the 
accumulator) is: set N\ and subtract it twice from the accumulator. The 
first subtraction gives zero, and checks the setting: the second giv^ JV, 
together with some irrelevant 9*s in the extra digit-places of the left 
half of the accumulator. (Try this.) 


THE DIFFERENCE TABLE 
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1.7. 

A dififereoce table, for a function at equal intervals. 

Serial No, Value 1st Difference 2nd Difference 3rd Difference 

0 -10 

1 - 4 

2 0 

3 +2 

4 +2 

Above is an example of a Miflferenoe table* for a series of values. When¬ 
ever such values are derived by putting successive integral values for n 
in an algebraic (polynomial) formula, the differences become equal all 
the way down one of the columns, and zero in the next, as seen here. 

Set out a difference table for nos. 0,1,8,27,64,125,216. In the column 
in which the differences become equal, put an extra such difference at 
the foot. Use it to construct an eighth number of the sequence given. An 
example of this procedure is shown in the note for Ex, 4.4. If you can 
also fit such an algebraic formula /(n) to the given values and calculate 
/(7), then you will have checked your work. Otherwise, the check for 
differences is to add back, e.g. in this table the first differences are 
checked by starting at the top of the ‘values* column, and then adding 
successively 6,4,2,0, getting at each stage the next number down in the 
value column. (It may save time in simple cases to subtract mentally to 
make the table, and only use the machine to check it.) 

1.8. Construct a cubic polynomial an® -f + cn + d in which a, 6, r, d 
are chosen integers (positive or negative). Calculate the values of f (n) for 
values of n from zero to 5 i.e. tabulate f(ii) for n=0(1)5; the bracket 
denotes the interval between successive values of n. Form the difference 
table and verify that the third differences are equal. Show algebraically 
that such equal differences must have the value 3 x 2 x a = 6a. 

[Solution note. 
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CHAPTER 2 


MULTIPLICATION 


Learn how 

(i) to multiply a given number, e^g. 37-48, by an integer such as 297; 
operating first with the first figure (2, here). 

(ii) how to do the same more shortly as 300 - 3, checking the process 
by the figures in the multiplier register. 

2.1. Multiply (i) 57-88 by 967 

and (ii) 9-67 by 5788 

Do each (a) by giving the full number of turns indicated by the 
multiplying figures (e.g* % 6, 7) and (6) by using backward turns to 
reduce the total. 

Note: Calculation can be done in as few as seven turns. 

Consider how you would have put the decimal point markers if the 
calculation (i) had been 578-8 by 96-7. Markers should appear (a) in the 
setting register, where the figures 5788 are set up; (b) m the multiplier 
register, where 967 appears after the calculation is done; and {c) on the 
accumulator for the result. 

Notice that the decimal point markers do not in any way affect the 
operation of the machine, but only the interpretation of input and output. 

For 2.2. If your machine allows it, learn how to transfer a result from 
the accumulator back to the setting register. 

2.2. Multiply 4*83 by 5-54 and the result by 22-3. Check by carrying out 
the multiplication in a different order. (If you have a back-transfer from 
the accumulator, use it to set your first product, ready to form the triple 
product.) 

2.3. Evaluate 3-79x5-4 + 12-6, first in this order; and secondly by 
adding 12-6 (in correct position) into the accumulator first. 

2.4. Evaluate 11 x 1-3 + 35 x 1-475 +24x 2-67 first with the integers as 
multipliers, noting the final figure shown in the multiplier register; and 
second with the decimal numbers as multipliers. 

Notice that although the latter method takes longer it gives a result 
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ACCUMULATION 5 

which might often be required viz.: the sum of the decimal numbers, as 
in the following example: 

Determine the total cost and weight of a freight consignment of 1-3 
tons at Sll per ton, 1475 tons at S3S per ton, and 2-67 tons at S24 per 
ton. 

2.5. The heights of 100 boys are measured to the nearest inch and are 
found to fall into these intervals: 


Inches 50-52 53-55 56-58 59-61 62-64 65-67 68-70 

No. of boys 3 12 19 30 21 11 4 


Ascribe to each boy, for convenience, the mid-height of the interval in 
which he comes; and work out the estimated total height 51 x 3 + 54 x 12 
and so on. 

Hence find an estimate of the average height. 

(Notice the check the multiplier register gives, that the numbers of 
boys have been correctly included in the calculation.) 
l^he result is checked by the calculation of Ex. 2.6. 

2.6. The data are as in Ex. 2.5, but the width of the interval (3') is taken 
as a unit of height throughout the work, and the heights are measured 
not from the floor but from 60* level as a standard; heights above being 
counted +veand below -ve, thus: 

Height in units -3 -2 -1 0 +1 +2 +3 

No. of boys _ 3 12 19 30 21 11 4 

Evaluate the total of heights counted in this way, and the average. 
Total is 3 X ( - 3) + 12 X ( - 2) + etc. 

(Notice that since we cannot set negative numbers but we can 
multiply by them, the multiplier register reading will not this time give 
the number of boys. Will its value help us in any way to check that we 
have not omitted data or Included any twice by mistake?) 

Remembering how to interpret the final average in terms of units of 
3*, above the 60* level, verily that your result checks with that of Ex. 2.5. 

2.7. (a) {b) 

Evaluate 2-431 x 4-7 - 0-09 x 10 8 

(i) by forming (a) in the accumulator first (but using either factor 
as multiplier according to convenience), 
and then forming(6) negatively; 

(ii) by forming (6) negatively first and then adding (a). 

& 


M.K.M, 
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2.8* {a) Evaluate 3-9 x 4*88 and 6*2 x 4-88 in a single operation — by 
setting the first numbers far enough apart on the setting levers, 

(Check; Verify that the sum of the results is 10-1 x 4*88.) 

If your machine has a locking device for part of the accumulatorj 
learn how to use it, 

(6) Repeat Ex* 2*8 getting one result where you can lock it* Do so; 

then set the other, shift so that the two are in alignment, and accumulate* 
If your machine has back-transfer from the accumulator to the set¬ 
ting-register, proceed as follows: 

{c) Repeat Ex, 2,8 (b), but transfer back the one result instead of re¬ 
setting, (The other result, being locked, will not transfer.) 

2,9, A number of parallel forces (or weights) etc. act at points 

on line Ojt, drawn perpendicular to them, at points distant etc. 

from O, The sum (M) of the moments about O is 2 (i*e* 

WiX^ + IFgXg + etc. 

and the sum (P) of the forces is i-c- 

W^g + etc, 

Find in a single series of operations the values of M and P and calculate 
5=M/P, 

Value of W 4 3 5 2 6 (kg.) 

Value of X 8*9 9*4 9-7 10*4 10*8 (metres) 

Repeat the calculation, where x is measured from the 10 metre point, 
so that the values of x read: 

-M, --0*6, ^0*3, +04, +0*8 

Verify that the position of the centre calculated in this way is the same 
as betbre. 


CHAPTER 3 


CONVERSION OF UNITS, INTEREST. ERRORS 

This chapter is an extension and continuation of the previous one* 
What k new is arithmetical rather than mechanical, 

3.1, (a) Find 173*24 x 7*9, Note, but do not clear. 

Now find 173*24 x 8*2 simply by making 3 extra turns. 

Note, but do not clear. 

Check by setting the hrst result and subtracting from the second 
(which stands in the accumulator), 

(6) Use the same method {continuation) to find 381*75 x 0*14 and also 
381-75 X1-14 and check similarly. 

Construct a commercial percentage problem to which these two 
results are the answers. You may prefer to set the problem in a country 
with decimal coinage. 

3.2, Write down a series of small sums of money (in shillings and 
pence). Set them in positions 5, 4 (shillings) and 2, 1 (pence) and add 
them successively into the accumulator. Rectify at intervals by adding 
the number 988, which subtracts 12 from the ‘pence’ and adds 1 to the 
shillings; and continue until the result is obtained. 

Devise a similar method of rectifying shillings which have exceeded 
020 (but not 999) in the section 6-5*4 of the register. 

Add a series of sums of money, in £.s.d. without fractions of a penny, 
and check by adding in reverse order. 

Devise a similar procedure for dealing with stones, pounds and ounces, 
and test it, 

3.3, Write down a series of lengths in terms of feet and inches* You are 
required to tabulate them in column, with metric equivalent in cm. 
(2*54 cm, = 1 inch) alongside. 

The centimetric column should have its total at the foot* 

Method, Convert all distances to inches and decimals mentally and 
note them in column. Multiply each in turn *by’ 2*54 — which number 
will be set on the levers? — and note the result. Clear only ibe ac¬ 
cumulator and setting register after each of these operations* What 
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CONVEHSION OF UNITS 


number will then appear in the multiplier register? Use it to get the 
final result required. 

In order to check, would it be enough to verify that this final answer 
is the sum of the separate answers? (You will have to consider whether 
errors which arise are due to you or to the machine.) If your answer is 
** no”, devise and carry out a satisfactory check, 

3.4. A journey consists of a scries of stages, each given in kilometres. It 
is required to set out, in the form shown by the similar example below', 
both the separate mileages and the running total in miles, thus: 



Distance 

Running total 

in km 

in miles 

in miles 

20 

12-40 

12-40 

30 

18-60 

31-00 

11 

6-82 

37-82 


Data for your example 20*8,33*4,17*4,9*3,18*5 km. 

Conversion factor 1 km, =0*62 miles. 

Method 1. As in Ex. 3.3. Convert the stages first to form col. 2, and 
form the running totals afterwards. Notice that you have room to set 
two of the data so as to multiply both in one operation. 

Check (final total) by taking data in reverse order or by addition of 
data in kms. and conversion of this total. 

Method 2. (Available if left-hand side of accumulator can be locked.) 
Clear the accumulator and lock the half. Set 208 twice on the levers or 
keyboard, so that it will feed into the accumulator at the right*-hand end 
of each half, e.g. at 00208/00208, 

Multiply by 0*62. This gives the mileage twice. When it has been 
noted, clear all registers (except the locked half which will not clear and 
will give the running total). 

Repeat with 00334/00334: the right-hand answer only need be noted. 
Continue to the end. Check as in Method 1, 

3,5. If Ex, 3.4 has not been done by the method shown, find with the 
same data in km, (20*8,33*4,17*4,9*3,18*5) the total distance in miles, 
taking 1 km.=0*62 miles. 

Find the greatest possible error which can have arisen in the result 
(i) because 1 km, has been taken as 0*62 instead of 0*621 miles, and (ii) 
because the data are measurements known only to the nearest tenth of 
1 km, [Solution note. 


INTEREST 
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Note these abbreviations 

2D for two decimal places, 

4S for four significant figures. 

A note on the rounding of figures (required for Ex. 3.6), 

The value 0*128, rounded to 2D, is written 0*13, except when it 
represents a sum of money; c.g.: a 10% service charge on a bill for 
T28 Francs, (Le, 1 Franc 28 centimes) is 12 centimes. The last figure is 
cut off, whatever it is. Since no correction is made, the process is quite 
mechanical: the percentage could be added by transfer-and-shift 
‘without looking*. 

In ordinary numerical work, however, the 2 would be made 3 by the 
operator, (Brunsviga ‘20* allow'S this, direct to the accumulator, by the 
thumb-wbeek.) 

3.6. A sum, deposited in a certain bank, yields 4% interest after one 
year, and this if not drawn out it is added to the original sum. Interest at 
the end of a further year w'ill be counted on the total sum credited to the 
customer at the beginning of that year, and so on (compound interest). 

The sum of SI247 is deposited. List the interest credited at the end of 
the year, and the amount; and so on (no money being drawn out) for a 
second and for further years up to ten. Set out as shown below and check 
by adding the total of interest payments to the original sum (or in detail, 
as for a difference table). This is the layout: 

Date Amount ,5 Interest $ 

Commencing date 1247*00 — 

End of 1st year 1296*88 49*88 

Notice that the result of the first stage of the calculation is the starting 
point of the next. If you are using transfer from accumulator to setting 
register, make sure to have the decimal points in the same position before 
transfer, 

3-7. Treating the table of Ex, 3.6 as the beginning of a difference table 
for a ftinction f{n)j the amount after n years, would you expect some 
higher order of differences to be constant — apart from the slight effect 
of rounding off figures? 

You could consider this problem algebraically, writing 

f(n) = P. R" where R = 1*04 in this case. 

If the figures of Ex, 3.6 give an inconclusive result, form a table for 
P = a large power of 10, and R = 1*1, so that rounding does not obscure 
the result. 

Contrast Ex, 1,8. Can you generalise the result given in that example? 
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ACCEPTANCE 


For 3*8, To determine the greatest possible error in certain functions 
when the arguments are subject to errors up to a given magnitude. 
Argument is a better name than variubk m this context, 

(i) If a function of x is tabulated, the result required can be read off ; 
c*g. if as = 304“, rounded to 3S, (i.e. between 304 ± 0*05) then the 
difference for 3 minutes gives the range of sin x to be 

0*5000 ± 0*0008 


and this states the acceptance of the value 0*5000, 


(ii) For simple cases the result can be written down. 
If = 4*3 ± 0*1 and y = +9-7 ± 0*2, 
then 

x^y= -5*4±0*3 


(iii) For other cases an algebraic treatment can be used e.g, for the pro¬ 
duct suppose the true values to be x-iAx^ y-^Ay\ then 
neglecting the term Ax , Ay the error is seen to h^x,Ay-\-y. Ax 


(iv) Failing any other method it is possible to vary the values of the 
data, through their possible ranges, and to see directly what is 
the effect on the result. This is in fact what we have done in (i). 
Note that in all practical contexts a numerical result, though fully 
checked, is of little value unless its acceptance is known. 


3.8, Given x *= 4*52, y — 0*73 (rounded values) calculate the values and 
acceptances of 

(i) tan y, {y in radians) 

(ii) x+y 

(iii) xy 

(iv) 

Check (iii), (iv) by direct calculation of an extreme value* 

[Solution note* 


FLOW DIAGRAMS 


11 


A note on Flow Diagrams 

In many instances we find that a specified calculating procedure is 
carried out again and again, each time with different data, until some 
prearranged condition is satisfied. 

A clear way of setting out such a calculating programme is by 
du^ram. We show such a diagram below, devised for the solution to 
Ex* 2.5 on p. 5* 

Instructions are contained in rectangles; and each rounded enclosure 
contains a yes-or-no question* 

In this example we use abbreviations: 

A for accumulator, C for cycle register. 



(Result is in C) 


In later examples of flow diagrams the machine-details may be 
omitted, instructions being given in arithmetical or algebraic form. 


FROM THIS POINT, ALL PRACTICE EXERCISES 
WILL BE FOUND ON RIGHT-HAND PAGES* 
NOTES FOR EXERCISES ARE ON LEFT-HAND 
PAGES 
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NOTES FOR EXERCISES 


For 4.2. (Example is on right-hand page). If you have a lock, do the 
example thus: set the integer in two positions before squaring* then one 
of the squares is cleared and the other is accumulated. 

For 4,4. 


n 

f(n)=n^ 




A*£ 

0 

0 


— 





1 




I 

1 


14 


— 



15 


36 


2 

16 


5{}^ 


24 



65 


60 


3 

81 


110 


24 



175 


84 


4 

256 


194 

(assumed) 



369 


108'^ 


5 

6 

625 

1296^^ 

67r^ 

302"^ 




The arrows show how, when a constant value for a certain order of 
diflFerences is regarded as established, the table can be extended: but 
periodically a direct check must be made, e.g. here first at 10*. 

For 4.4. (ij) The coefficients in these equations are those of the 

Bmomial Expansion 

(1 + 2 )" = 1 + + + ... 

For integral exponent Ji, they can be built up for successive values of n 
by Pascal's Triangle 1 1 

1 2 1 
13 3 1 

where each figure is the sum of the ones above and ‘to the North-West', 
giving next: 1 4 6 4 1 

(b) For non-integral exponent the coefficients must either be cal¬ 
culated as 

P P(p-n P(p-^)(p-2) 

1 2! 31 

etc. or obtained from tables of interpolation coefficients. 


CHAPTER 4 


SEQUENCES: EVALUATION. SUMMATION 
AND DIFFERENCES FOR POLYNOMIALS 


4.1. A sequence of numbers can be evaluated in turn, 

given the value of Sq and a rule connecting successive values. In this case 
the rule is S,, S„_i(0*8) +1 and So=0, 

Find the value of to 5D, 

Verify by finding an algebraic formula for (or failing that, for 5^) 
in as concise a form as possible. 

(It may be convenient to write r for 0*8 and to substitute its value at 
the end). Hence check your result, [Solution note. 

4.2. Sum the squares of the first ten natural numbers (and record the 
individual values if it is possible to do this in one operation: sec note). 
Do not char tlie muitipUer register between successive squarings. 

What is the number accumulated in this register? 

Verify this value and the sum of squares by using the algebraic 
formulae for 

2(r) and forn=10. 

r« 1 r— 1 

4.3. Find by direct accumulation the value of 

302 - 292 + 28^ _ 27^ +... + 22^ - 21^ 

Verify that the value is 10(30 - 5 - 21)/2. Why? [Solution note. 


4.4. Set out in column as shown the values of/(#i) = n* for ir ^0(1)10 and 
differences. The top differences are regarded as linked to the first value 
of n, as shown, ♦ and are written as ^/(O), jd®/(0) etc., e.g. -dy(0) = 36. 

(a) Verify that /(2) =/(0) + 2J/(0) + JY(0) 

/(3) =/(0) + 3d/(0) + 3 jy(0) + AJ(0) 
which are examples of the Gregory-Newton formula 

f(n) =f(0) + (j)^f(O) + ... + (“) J'f{0) +... . 

Clearly this would not he used to find/(«) for integral values, but 
it is used for/(^) where p is not integral. 

(b) Use the formula /C/»)=/(0) + ^^^J/[0) + ..., to evaluate^i). 


First list the values of 



Check that the result is the fourth power of 
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NOTES FOR EXEECISH5 


For 4.6 (a) represents *?(?» +1); + !)(« +2) and so on, 

while =if. (This example is for the more advanced student,) 

(b) We can use the Gregory-Ncwton formula, with n integral, for 
summation of a polynomial sequence* 

If we write 2 /(O ^ fi^) = ^ S{k), 

F —0 

We regard the/column as the first difference column of the function S. 
Then 5(10) = 5(0) + + (^2 the 

terms so long as the differences are non- 2 ero. The term S(0), as the sum 
of zero terms of the/sequence, is zero* [See Finite Differences, p. 65.] 

For 4J. The p esmt value of a payment P a year he* ^ is 

P/14W = (0%15)R 

(Treat 0*9615 here as if it were exact.) 

The present value of P to be paid n years hence is (0*9615)'* P, 

For 4.9. If a table of values of a quartic polynomial were required e,g. 
for x=0(0'l)l, it would only be necessary to calculate until the 4th 
difference is established, (see Ex, 4*4.) It is {^\)a^ for unit interval, and 
more generally (4!)/ioA^ for interval h. The h* factor could be verified by 
rewriting the Uble of Ex. 4.4 for 0-1, 0*2 etc*, which gi ves the same first 
column figures multiplied by 10“*, In practice the difference figures 
would be written without decimal points, the last figure in every entry 
being understood to be in the 4th decimal place. 

For 4.10. Detection of errors by dififerencing. 

If a table of values is known or believed to be those of a polynomial (or 
a function which can be so represented in that range) errors can often be 
detected by forming a difference table. 

An error ( + e) in a single value will affect the differences according to 
a fan-pattern thus: 


Error in f 

in At 

in AH 

in AH 




■he 








-3e 

+ e 

-e 

— 

+ 3e 




-e 


FUNCTION VALUES 15 

4,5. Sum the cubes from 11® to 20®. Check by using the formula 

2 (^^) ^ (n + 1 )*/4 for n =20 and for n = 10. 

1 

4.6* (a) Establish in general (or for r « 2,3,4) that if ff(n) is then 

4/‘r(«)=/r-l(n + l)- 

Use this result to build up in turn the values of / 3 (n), fjn) and f^(n) 
fromn=0 to 10- 

Check directly for u = 10. [Solution note, 

10 

(b) Verify the equation of note 4,6 (b) for 2/i(^0‘ 

n“l 

4.7, Find, to the nearest cent, the present value, reckoning interest at 
4% compound, of a series of 5 annual payments of $ 550, the first 
payment being a year hence. Before deciding your method, refer to 
Ex* 4.1. 

Check by formula he. 25P{1 -(0-9615)^} 

4.8. Evaluate to 3D /(x) = afpe* + where = TOO, 

Ej = 0*87, a^ - - 0 61, Eg T20, a^ =* 2*00 and at = 1-IS. 

Method L 

List the powers to 5D. (When has been found, leave it set to find 
jc®, X*). Evaluate term by term. 

Method IL 

Find W|=ao.r + a|, then W 2 = WiX + % and so on. Finally is the re¬ 
quired result. 

Find the greatest possible errors in the result in cases (i) x supposed 
exact and the a*s ‘rounded^ e.g* 1 *00 means some value between 0*995 
and 1*005; and (ii) the a's are exact, and x subject to an error ±0*01* 
(See Ex. 3.8.) [Solution note. 

4.9, Construct a table of f(x) = 2‘50x*-f 10‘7x^ -0*89 for x =^0(0*1)! 
using method II for x = 0(0*l)0*k 

Build up thereafter by differences; and check/(l) directly. 

4.10. The following are gi ven as values of a function at equal intervals of 
the variable: 

0*01, 0*01, 0*13, 0*73, 2*41, 6*10, 12*61, 23*53, 40*33, 64*81, 99*01, 

145-21. 

One value is in error; discover and correct it. [Solution note. 
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NOTES FOR EXERCISES 


Learn how to divide by ‘teardown divbion’j which is the long 
division method of paper arithmetic: good examples are 

(i) 1 'OOGO by 1, What is the result to 4D ? 

(ii) 12-341 by 0-25 exactly. 

Note particularly the locations of the data and the answers — quotient 
and remainder — in this process. 

(a) The dividend (N) is required to be in the accumulator at the kfukand 
end of iL It is set, and transferred to the accumulator by a forward turn: 
the figure 1 which the turn causes to be shown in the multiplier register 
must be cleared immediately^ and N itself is cleared from the set. 

(^) The carriage is now shifted fully to the right> and the divisor (D) is 
set in the correct position for subtracting from the highest figures of N, 
[The object of (d[) and (^) is to ensure that the first quotient figure will 
appear at the left-hand end of the cycle register, and thus enable its full 
capacity to be used.] 

(c) Division can now begin. 

The whole process is shown in a flow diagram: 


Carry oui (c?) and (i) 


Crank backwards once 
Q Beil rings ? ) 


-NO 


YES 

I_ 




Crank forward once. 

STOP 




Push carriage left one position 


is in cycle register) 


-NO 


YES- 


CHAPTER S 


DIVISION* SUMMATION OF POWER SERIES 


5.L (( 2 ) Divide the answer to Ex. 2.2 x ^ x r) by each of its 3 factors: 
there will be no remainders. In each case check the decimal point of 
your answer by inspection of the product of the other two factors, i.e, 
test (pqr)lr by inspection of^^, 

(b) Find the reciprocals, to 4S, of 1-249, 0-G8334, 750-2. Check the 
figures of your answers by reciprocal tables, and then multiply each 
number by its 4S reciprocal. 

Can you see why the discrepancies differ in magnitude? 

[Solution note. 


5*2. Calculate 


2 47x39-2 
43-12 


(a) in order, as you would read it; and 

(A) by a method suitable if a series of numbers were to be reduced in 
the same ratio 39-2/43-12, 

What is the greatest possible error in your result due to rounding errors 
in the data 2-47 and 39-2? 

(It can be seen that the greatest possible error due to 43-12 is smaller, 
and you are invited to neglect it here.) [Solution note. 


5*3. (i) An alloy is made from 10-3 kg* of copper, 27*4 kg. of nickel and 
0-5 kg. of tin. What is the percentage of each metal in the alloy? 

Check that the separate figures total to 100, to the accuracy which is to 
be expected. 

(ii) With the data of Ex. 3.8 find the acceptance of ;^ = xfy. Verify that 
it agrees with the formula 

(dz)lz ^ {Ax)lx + (Ay)ly^ all numerators positive. 

Show that the same formula applies to Ex. 3*8 (iii). [Solution note. 
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18 NOTES FOR EXERCISES 

For 5,4. Leam how to do 'build up' division, taking as examples again 


I 


I 


(i) 1-000 by 7*0 

(ii) 12-341 by 0-25 


This method is preferred when the 'turns' can be used at the same 
time to multiply another number, e.g. in (ii) the number 2-43 could also 
be set (elsewhere on keyboard) and w^ould give a result in the ac¬ 
cumulator equal to: 


2*43 X 


12-341 

0-25 


For 5.5 (i) If build-up division is used, notice carefully in which 
register the quotient is shown. In this problem the quotients are being 
accumulated, and this register must not be cleared after the first division. 

(ii) Notice that teardown division can be used cumulatively if you 
look after the 1 which appears in the multiple register when a dividend 
is set in the accumulator. This 1 cannot be cleared after the second 
numerator is set. 

It could be removed by making a back-turn before the figure is set. 


For 5.7- (i) A possible method is to tabulate n and l/(wl) for = 1 (1)15, 
to as many figures as the machine allows. 

(ii) Could you shorten the amount of paper-work by accumulating 
reciprocals as they are formed? See Ex. 5.5 (i). 

A flow diagram forjthis work; 


Take sum ^2-5; term Tr = V6; r = 3 



1 . 


DIVISION PROBLEMS 
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5.4. Repeat part of Ex. 5.3, i.e. find the percentage of copper by 
'multiplying up* the total weight to make the weight of copper and at the 
same time ‘multiply* by 100 to give the percentage. 

(As stated, this setting of 100 is unnecessary since result figures are 
given in the multiplier register; but if the point has been followed you 
can repeat the procedure to get in one move the answer to the following 
problem: 

In 4.39 kg. of the alloy^ how much nickel would there be, by weight? 

(38). 

Check your result by using your calculated % of nickel from Ex, 5.3, 

5.5. Evaluate (if possible without getting both quotients separately and 
adding) 

1*74 2M 

6-231 ^ 1-75 

(Consider W'hich method of division to use) 

Check by repeating, forming the second fraction first, 

5.6. In one year a country exports 57,100 tons of a commodity for 
£238,400; and in the following year 51,500 tons for £253,100. Find the 
% increase or decrease in each case and write them to the correct num¬ 
ber of significant figures. Find also the % increase in the price, and show 
that the results are consistent. 

5.7. Devise a method of evaluating e correct to 5D (working to 7D) by 
means of the series 


+1/(1 l) + l/(2!) + l/(3 [) + .... 

Record your method in detail, and check your answer with the known 
value of e. 

It is also checked by Ex. 5.8. 

A self-contained check for this exercise, if method (i) is used, is to 
evaluate also 


^-^ = 1-1/(11) + l/(21)-l/(3!) etc., 

for which all the terms are known. Then check 1 to the expected 

accuracy. 
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For 5.8. 

Since x is the reciprocal of an integer, your method of Ex. 5.7 could 
readily be modified; but you should consider also how to programme if 
;c=0'51 say, 

(If you list powers, N.B. method I of Ex. 4.8; can be set to get 
values x^ to x^, but cross-check by e.g. x 

For 5.9. The result, log 10 == 2*3026, could be checked by evaluating e 
to the power 2*3026, A method would be 

(i) to find directly. 

(ii) to find by expansion and cube it 

(iii) and by a (very short) expansion 
and (iv) multiply all together. 

How many figures to take in each part, and how to assess the final 
figure ( = 10 nearly) are matters of nice judgment. 


For 5.11. Notice that owing to the nature of the number a: = ^ you will 
probably prefer to find powers by multiplication for the first expansion, 
and division for the second. 


POWER SERIES 
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5.8. Modify your method of Ex. 5.7 to find Je by the series 

= i +x+x®/(21) 4'X®/(3!) +, with x=-^. 

Get the result correct to 6D. By squaring this result, to as many 
places as your machine allows, check against the value of e. 

(Record the values of all terms of the series for use in Ex. 5.10.) 

5.9. Using the scries (valid if -1 <x< 1) 

=2{x + *3/3 +x^is + x3/7 + 

evaluate (i) log^ 2 (ii) log^ (5/4) and hence find log^ 10. 

(Consider carefully whether to get the later odd powers by dividing by 
Ijx^ or multiplying by air®.) 

5.10. Using the values of the terms found in Ex, 5,8 determine the sine 
and cosine of an angle equal to i radian from 

Sin X=X “ x®/(3 1) + x®/(51) - etc. to 5D. 

Cos X = 1 ^ x®/(21) + x^/(4[) - etc. 

Square and add, to check these results. 

Also convert ^ radian to degrees, using 7r = 3*14159 and compare with 
table values of sine and cosine. 

5.11. Evaluate tan-V(l/S) and tan-^(1/239) in radians to 8D, using the 
series, valid for -1 <af<l, 

tan^i X=x -x®/3 + x^/5 - etc. 

Hence find tt, usmgK = 16 tan-i(l/5) -4 tan^‘(l/239). 

Describe your method, and state with reasons how many figures in 
your result you claim to be accurate (or better still, the uncertainty in 
the last figures obtained from your calculation). 


c 


M.N.M. 
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A NOTE ON CHECKS 


Great emphaisis has been laid, in all the foregoing work, on the 
checking of results. This will have taught the reader those ways in 
which he is most liable to error, and (in places where the check show-s an 
inevitable discrepancy) the limitations inseparable from the results at 
the best of times; but it is clear that there is a heavy duplication of 
labour in, for example, the checking of an isolated result like that of 
Note 5.9- 

Two principal ways in w^hich this work can be avoided are: 

(1) Values of a function/(.v), computed for regular intervals of the 
variable, e.g. sin at 10 “ intervals, can often be made to check each 
other by examining the successive differences. Functions which can be 
represented by polynomials over the range of the tabic will be found to 
have a difference column which is constant as nearly as the rounding off 
of the values of /(at) will allow. 

The method of detecting a singk error in a table has been set out in 
Note 4.10. 

(2) For isolated values it may sometimes be possible to check by a 
quick process e.g. the root of an equation by substitution, but it may be 
hard to judge whether the inevitable discrepancy, the *residuar, is small 
enough — see later Note 8 . 6 . 

Wherever possible an isolated value is found by an iterative process in 
which an approximation is successively refined again and again until it 
emerges unchanged from the process — to the required degree of 
accuracy. The next chapter is devoted to tht^e methods. 

Square root. 



CHAPTER 6 


ITERATIVE METHODS, INCLUDING 
STARTING PROCEDURE 

Square root of a number N, 

To find the square root of a number N to any required degree of 
accuracy, starting with a reasonable first approximation Xq. Calculate in 
turn the numbers Xi=i(x 0 + iV/xo); X 3 =i(xi + JV/xJ and so on, the 
general relation being =i(^n ^l^n)* 

After a few iterations the successive values of x are found to differ by 
extremely small (and diminishing) amounts. When successive results 
agree to the degree of accuracy required, the answer has been obtained; 
e.g. if and Xg, rounded off to SS, are both equal to T4142, then this is 
the required square root to 5S. 

The method can be illustrated graphically by the sketch of the graph 
of y^Njx on the opposite page. 

The points Pq (xq, Njx^) and QoiNjxQ^ x^) both lie on the graph and 
are imag^ of each other in the line x =y, on which their nud-point Mg 
lies. The required root x=JV/x is the x-coordinate of the point F, and is 
equal to OR in the figure. 

It will be seen that the x-coordinate of Mg is the mean of those of Pg 
and i.e. ^(xg) +i(iV/xo), which is x^. We therefore drop a perpen¬ 
dicular from Mg to the x-axis, to cut the curve in Pi{xj, iV/xj whose 
image is Qi^ and the midpoint Mi ofPiQi gives in the same way P^; and 
so on. 

It is clear from the figure that the successive perpendiculars tend 
towards the position FP, and the x-values to the required root. 

Verify, on the same figure, that the values also converge to JN if Xg 
is smaller than the true value. 

(This analysis shorn the values x^, etc., as steadily decreasing; but, 
as in all digital processes, rounding may effect this unpredictably.) 

6.1 Use the method described to evaluate JIO to 4D starting with 
Xg « 3 -0. Note your programme of operations in detail. [Solution note. 
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NOTES FOR EXERCISES 


Note for 6,2, Newton *s Method for finding the root of an equation . 

This is an iterative process of very wide application which enables us 
to obtain to any required degree of accuracy the value of a root of the 
equation f{x)=0, given a sufficiently close approximation to this 
root. 

We write * 1=*0 * 


where f'{x) is the first derivative of/(*); and in general 
=*n - * ^®wto“ ’* Formula 



The point Pj is [aco, /(*o)] and the gradient at this point is 

T,P,ITM'{x,), 

Thus OT,=x,~jj^^=x, 

and a similar construction gives 0 X 2 =and so on. 

Notice however that if we start at 00, the turning point, the 

next result is worse. 

It will be seen that the vanishing of f*{x) at a point near Pq will demand 
more accuracy in the starting value It is such considerations, as the 
graph above will make clear, which indicate what is involved in the 
words ‘sufficiently close*. 

N.B. It must not be supposed that Newton*8 formula fails for a multiple 
root: the difficulty occurs for simple roots which are nearly equal, or for 
a root with a turning point very near it. 
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PROGRAMME FOR NEWTON*S METHOD 


6.2 (a) Notice that the formula given by Newton may require to be 
simplified before it is actually used for calculation, e.g. to obtain a square 
root we may write/(;r) - AT, which gives 

® 2xo 

Simplified algebraically this gives the formula used in Ex. 6,1. 


6.2 (b) Obtain an alternative formula for the square root, by writing 
f{x) = l/x® -1 /iV=0, in the form: 



A simple machine programme would be as follows: 

First ^culate and note the reciprocal of 2JV: clear. 

iteration (i) form the product ;c,j(l/2N)=P; 

(ii) set P and clear accumulator (manually or by direct 
transfer if available | 

(ill) multiplyPby {-x^:do not clear; 

(iv) add 1-5, to get value of bracket B; 

(v) set B and clear accumulator as in (ii); 

(vi) multiply P by to give x^^^ = X: note this value; 

(vii) set X and clear accumulator. 

Thereafter the process is repeated as often as is required. 

Use this programme to solve Ex. 6,1. 

Compare the two procedures for (a) convenience, and (5) rapidity of 
convergence. 

It will be appreciated that a machine can be made to carry out the 
whole process unaided if it has file following facilities: 

(i) for storing numbers which are required to be noted in the 
method above; 

(ii) for feeding in or tranrferrir^ numbers, as required for process¬ 
ing; 

(iii) for testing whether x„ and are or are not equal within the 
required limits, and giving an appropriate YES or NO signal; 

and (iv) a control nmV, to ensure that the prearranged programme is 
carried out, and repeated so long as the signal is NO, 

See also p. 58: Programming and the Computer. 
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NOTES FOH ESEHC13ES 


For 6.5. 
Formiilae are: 


^n+l 



and 



For 6.6. 


Iteration by direct formula for equation x=F(x) 

This works, writing =F(jc,») provided that the gradient F'{x) in 
the region used, is numerically<l. This can be seen by drawing figures 
as below for various cases. Sketchy=x andy =F(:r). 



The graph shows how the successive approximations are related to 
one another. It is clear that this is not an accelerating convergence like 
Newton's, which frequently doubles ike number of significant figures at 
each iteration. By contrast, this method gives steps which diminish in 
nearly constant ratio. 


ITERATION AND EXPANSION COMPARED 
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6.3. Derive from/(:K) =• l/« - iV a formula and a programme for finding 
a reciprocal, and use it to find 1/31 to 6S, starting %vith Xq = 0'03. Notice 
that division is not used. 

6.4. Starting with jeo = the value of square root of It given by 4-figure 
tables, find and using either formula from page 25. State the best 
result obtainable from this work. 

6.5. Using either the method of Ex. 6.1 or of Ex. 6,2 to get a Newton 
formula for a root, find to 4S the cube root of 28. {x^ = 3). 

Without completing a second calculation, consider the advantages of 
the above method compared with the use of the Binomial expansion: 

28i/3=271/»(28/27)^/3^3 (1 + 1/27)V3 

-K'n(r,)4.G)(-D{A)’-} 

6.6. Find a positive root of the equation x = 1-2 sin x (* in radians) to 
4S by direct iteration, using the table of sin x for x in radians, and the 
formula 


Xti+i = l’2 sinx„ 

If you have only a table of sin x for x in degrees, you could* use 
=sin x„ to get a root in degrees, and convert to radians. (Notice that 
4S is almost the same accuracy requirement in radians as in degrees. 
Why not exactly the same?) 

•It will illustrate the method equally well to solve the (different) 
equation « = 80 sin x to 4S where * is in degrees; but you should have 
access to a book of tables which includes functions of a radian argument. 
An excellent one is Attwood: Advanced S-figure Mathematical Tables 
(Macmillan). 














28 NOTES FOR EXERCISES 

For 6.7. To find a first approidmation to a root of f(x) 0. 

(I) Method of interpolation. (See footnote*) 

By trial we find values {X-[^AX) such that/(jp) has opposite signs 
for the two values. In the figure/(A')is -Pand/(^ + dA') is Q. 



x~X 


X~a x~(X^ 4 X) 


It will be seen that the chord intersects the x- 2 xis at a better approxi¬ 
mation {Xf) than either XotX-{- AX, 

The gradient of the chord = — 

Xf^X AX AX 

Xt -X=Pi^XIAf) = ~f{X)[AXIAn 
x,=X-f{X)[AXIAf] 

a fonn which will suggest the programme of computation. 

(H) Comparison of the diagrams for 6.2 and 6.7 shows that if x^ is 
too small, the value obtained by one Newton iteration from 
X or X + AX is too large; while a downward concavity 
reverses this situation. This suggests the rule: 

(a) From X and X-^AX calculate the interpolate value Xi as 
above; 

{b) from the better of X and X + AX calculate by one 
preliminary iteration; 

and (c) take Xq = J {x^ + x^) as the start for the main iteration. 

* Interpolation is the calculation of an intermediate function value. 
This example is of linear^ irwerse interpolation. 'Linear* explains itself, 
and 'inverse’ is clear if we consider the solution of say sin 0 — 0*6 = 0; we 
might estimate a value of 0 ^sin-^(*6) from 3 in-'(* 5 ) and sin”^(*7071). 


EVALUATION OF HOOTS 
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6.7. It will have been noticed that an iterative method is very good once 
it has been well started. A good start may save several iterations, 

A carefully drawn graph gives a good start, but judiciously calculated 
function values in the right region can be used without graphing by the 
method of interpolation (see note on left-hand page). 

Example: To find, to 4S, the real root between 3 and 4 of the 
equation x^-3x-24 = 0. 

We could take/(x) for which /(3)=-6,/(4)=28 

but with steadier gradient isg(jt:) = 5 ;® - 3 - 24/^c. 

Interpolate between 3 and 4 by method (I) of left-hand page and use 
this as starting value for Newton, with either / or g, (ii) Repeat the 
solution using method (II) to find a starting value. 

6 . 8 . Find to 4D a solution of e* =3x, given that: 

giMoo ^4-4817 
eLwoo =4-5267 
eX*62oo =4-5722 

Method J. Let x = 1*51 + h, and remember that = ct-sie* 

Method //. Work as in Ex. 6.2 (ii) for/(a!) = e® - 3jf. 

Here as in Method I, values of e® not given in the table must be found 
by expansion. 

The check is in any case by expansion. 

Notice that if the function / were known only by its table we should 
have to make assumptions about its behaviour which could only be 
tested by higher differences. See Ex. 4.4 ( 6 ); and for further reference: 
Interpolation and Allied Tables (H.M. Stationery Office). 

Method HI, See Ex. 6.6 and Note 6 . 6 , 

In the direct form this method fails to converge. 

The form Jfn+i (^^n) can be used, starting with Xo = l’5089 
derived from the value log^ (4-5267) = 1*5100. 

The reversed equation is, however, better solved by Newton’s method 
with 

/(a[:)=a:-log,(3a:). 

Use log*(fl + h) -= logjU + +etc. 





















CHAPTER 7 


STATISTICAL AND ALLIED CALCULATIONS 

Statistics is a branch of applied mathematics which is concerned with 
deriving from a set of numerical data —often a very large set — 
certain numerical measures which will indicate the principal features of 
the set, either considered alone or in comparison with another set* 

It is helpful sometimes to remember that the word statistics can be 
considered as a plural* A 'statistic' is a figure worked out from the data, 
c.g, from a census* It might for example be the proportion of males In 
the population of the country; or simpler still, the total population 
itself* 

A set of numerical data usually consists of the measurement of one 
characteristic of every individual member of the group: it could be for 
example the height (and possibly the weight as well) of every man in a 
certain group called up for military service, 

(1) The first statistic is the number of individuals in the group 
(symbol JV)* This should always be recorded, 

(2) The second statistic is the mean value of the measure which is 
being considered, e,g. the mean height of the men. If Xi is the height of 
the first man, X 2 of the second, and so on to for the last; then the 
mean value is 

( 1 ) 

The formula is sometimes used in this form, but usually two modifica¬ 
tions arise: 

(a) The measurements are made to only limited precision (c.g, i", in 
height) and for convenience we may wish to group results to an even 
wider interval. In either case a number of men belong to the same cla,ss 
(e*g, 23 have x = 28 have where 65 means a height of 65 

or above but lower than 66), The number 23 would be called the 
frequejicy of the value measured as 65* 

Then «= (l/Aot (/^,) (2) 

where is the value and /j the frequency for the first c/aii, and c is the 
number of classes. This ciculation is much shorter, but presupposes the 
classification of the data — a process which serves another important 
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purpose, viz,: to put the data in a form (pictorial or otherwise) in which 
their general nature can be exhibited to the eye. 

The shape of the data is often best shown by having about nine groups, 
e.g, for heights of men (i) under S' (ii) 5' but under 5'2' (iii) 5'2' but 
under 5'A* and so on to (ix) 6'2' and over. The mean obtained by group¬ 
ing in this way, with 10 terms in equation (2), would be slightly different 
from that given by using ail the possible groups, but would be clearer to 
see and easier to work with, 

(Another consideration affecting the mean is discussed on p, 32.) 

7,1, The following are the heights to the nearest inch, of a series of 
individuals. Classify them {a) into inch classes; and (ft) into classes 
50-52 inclusive, 53-55,56-58, etc. 

Values of height: 

55.57.60.64, 52, 60,68, 62, 57,58,55, 61, 62,58, 54, 69,68, 62, 50,59, 

66.57.61.65, 54, 61,60, 58, 54,63, 70, 59,54,57,53,64, 56,61, 63, 54, 

64.62.58, 65, 57,58, 61,66,66,59, 61, 56,60, 65,64, 60,67,64, 62, 65, 

59.61.59, 63, 62, 68, 58,61, 53, 64,59,62, 50,59, 67,61,61, 58,59,62, 

56.59.64.60.65, 57, 58.61, 60, 54,53,60, 55, 59,63,59,62, 57,64, 56* 

It will be necessary to use a systematic method of classifying, c.g, to 
give a column to each class and to enter each number in turn as a stroke 
in its column (as well as crossing the number off the list). It is con¬ 
venient to put each fifth stroke across the rest (thus, 44+4) so that the 
pattern of strokes at the end indicates the count clearly and accurately* 
The stroke total should agree with a direct count of the data (but this is 
no safeguard against entry in a wrong column)* 

After classifying, calculate Je (a) from the inch-classes, and (ft) from 
the 3-inch classes. Notice the small value of the discrepancy, (In case 
(6) the value of x taken for the class is to be the mid-value of the class* 
Sec Ex* 2,5*) 

Refer now to Ex. 2.6 for a simpler method of calculating by using a 
class-centre as origin and the class-interval as unit* Repeat this, taking 
values of x from 0 io 6 instead. 
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The calculated mean: dependence on the type of measurement* 
It should be noted that the calculated mean has to be understood with 
reference to the treatment of the original measurements. If the measure¬ 
ments were taken disregardifig all fractions of an inch the calculated figure 
is too low by an amount which can be taken as 4^5 but if they were to the 
nearest inch it can be accepted as it stands. 


For 7*3. 

The third statistic which we shall consider provides a measure of the 
‘spread’ of the data. Similar electric lamps made by two makers P and Q 
might prove, by a test of two equal sample groups of 100 (collected in 
some random way), to have the same mean Ufe; but if P’s output con¬ 
tained appreciable numbers (i) with less than one third and (ii) with over 
twice the average life, then this pattern would be less satisfactory because 
less dependable. 

It would be unreliable to measure the difference between the best and 
worst lamp, A possible method is to measure the difference between the 
25th and the 75th in order of merit in each group; but the method 
preferred is to find the standard deviation (a) of each group. 

We compute the deviation of each value x from the mean (j 5) ; sum the 
squares of these deviations and divide by iV, the total number of in¬ 
stances, to get the variance V. 

V (1/N)S (x “ 5)“ summed over all instances; 
but since there are /instances of each value at, we in fact compute 
V=(1 /N)S{^r ” 5)*} summed over the number of classes. 

Then o is defined to be 

This is the true standard deviation. 

It will be seen that the direct calculation involves all the separate 
deviations being measured from J?, a non-integral value. An indirect 
method will be seen later which avoids this complication* 
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7*2* 

You are given the following accident figures for a year in a certain 
district: 


n = 

0 

1 

2 

3 

4 

5 

6 

7 

8 

>8 

No. of days in 
which n fatalities 

occur 

30 

99 

IIS 

81 

28 

9 

3 

0 

1 

0 


Find (i) the total number of fatalities and (ii) the mean number of 
fatalities caused per day. 

Notice that the original data were a series of dated records each 
stating the number of fatalities caused on the day in qu^tion. Classified 
they appear as above, the frequency / being the number of times the 
number n appears in the records* («is used in place of x since it is of its 
nature integr^, though this not true of H=(l/A/')2 (f/r)-) 

r—0 

The summation is carried far enough to cover all cases: as always 
iV^=S/r over all classes (i,e, including/o). 


7*3, Referring to Ex. 2,5 and Ex, 2,6 (which is a restatement in more 
convenient form) we will show how to arrange a calculation of a 
directly. 

Herein+0 03 

Set out the values of x-£ for each class, with the frequency of the 
class; preferably in columns with the latter first: 


f 

x-x 

f(x-S) 

f{x-x)“ 


3 

-3-03 

-909 

27-54 


12 

-2-03 

-24-36 

49-45 

Complete the last 

19 

-1*03 

-19*57 


column: sum it 

30 

-003 

-0-90 


(B): Calculate 

21 

+0-97 

+20-37 


V=S/N 

11 

+ 1-97 

+ 21-67 


and a=VV 

4 

+ 2-97 

+ 11-88 


This is in the units 





used for measur- 

= 100»JV 


Sum=0*00 

S = 

ing XI convert to 



(check, not 


inches. 



usually exact) 
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NOTES FOR E3LERCLSES 


Standard deviation using an origin other than the mean. 

The same drill is carried out vrith values calculated from a more con¬ 
venient origin at a distance b from the true mean* The sum S of the 
squares so obtained is a!w£^^ larger than from the mean. The value of a 
is obtained from the equation 

a^ = SIN~b^ 

This result is exact, and is used in either of two ways, according to the 
chosen point for measuring the deviation of each value: 

(A) using the cetUral class-position as origin: this has the advantage 
that the correction is rather small — it may even be negligible 
to the accuracy required, 

or (B) using an end class-position: this has the advantage that all the 
numbers in the calculation are positive* 

Examples: (A) Ex, 7* 3 gives, taking zero as the origin 


f 

X 

fx 

fx* 

3 

-3 

-9 

27 (6--0-03) 

12 

-2 

-24 

48 

etc. 




Sums N «100 


+ 3 

195 = S 



=NS 



Then S/N = 195/100 -1*95 and = 1-95 - 0*0009, 

(B) Using the end class as zero, as in the Ex* 7,2* (Write x for n) 


f 

X 

£k 

fx* 

As X is already measured 

30 

0 

0 

0 

from the end we need 

99 

1 

99 

99 

not change it. 

115 

2 

230 

460 


81 

3 

243 

729 


28 

4 

112 

448 

False <i®=2133/N 

9 

5 

45 

225 

and x=75S/N=b 

3 

6 

18 

108 

Hence complete to find 


0 

7 0 

0 true o* 

1 

8 8 

64 

366=iV 

755 

2133 
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7 . 4 , 

Find the standard deviation of the following dassihed data, (They 
represent the weight distribution of 100 boys, grouped into 7 classes,) 
The value associated with each class is given the symbol y* 

y —3 -^2 — 1 0 +1 +2 - 1*3 

/ 3 11 21 28 25 10 2 

Obtain c by two methods, 

(i) with the origin as shown: and 

(ii) moving the origin to the left-hand class so that y has values from 
0 to 6: you should check on obtaining y in case (ii) that it is con¬ 
sistent with the position found for the mean in (i), 

7.5, Check your value of ct® for the accident figures on the left-hand 
page, by recalculating with x measured from an origin near its mean, so 
that X reads down the column from - 2 to + 6 i,e, method (A), 

7.6, The superiority of method (B) for machine-work is shown by the 
following programme for the accident figures. 

(i) Set the first value of x in position 6 and in position L (Positions 
are counted from the right in the setting register,) 

(ii) Make a number of forward turns equal to the first value of/, 

(iii) Clear the setting register only. 

/is now shown in the multiplier register, and in the accumulator 
fx and/:r® separately. 

(iv) Continue as above, with the next value of x and its square, and 
the corresponding /* These are therefore accumulating, and all 
the required sums are obtained together* 

Do this with the figures in (B) opposite and verify the column totals. 
Do the same with Ex, 7,3 taking 0 to 6, 











$6 NOTES FOH EXERCISES 

Far7J, 

Product-moment correlation. 


We have considered already the distribution of height and of weight, 
in suitable class-intervals, of a certain 100 boys. 

It would have been possible from the data to exhibit the distribution 
in respect of both attributes at once, as shown in the grid. Ex. 7.7 where 
/ (the figure in any location of the grid) shows the number of boy^ who 
have both the height and the weight classification corresponding to that 
location. 

We should expect the boys who are above average in one respect to be 
so also, more often than not, in the other; i.c. we should expect + values 
of X and y to go often together, and also - with -. Thus the sum of (ji^) 
products over all the individuals, if strongly positive, would indicate a 
positive correlation and could be used as a measure of this: we therefore 
calculate ^.nd (1/N)2(ficy) is called the product-moment 

about the origin i, e. about the centre of class (0, 0). 

This calculated result needs to be corrected because the origin used is 
not at the true mean: the corrected value is 

(1 (fxy) - ify (to the true mean) 

(Note that the correction term £y is not necessarily positive in this case.) 

This result is of more value if it is divided by a* and by the stan¬ 
dard deviations for height and weight (expressed still in terms of class- 
interval as a unit). 

The final result is called the product-moment correlation co¬ 
efficient of the data. Calculated as described, it is bound to lie between 
+1 (when all the results lie in the NW-SE diagonal locations) and -1 
(for the other diagonal). 
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7.7. 

HEIGHTS (x) in intervals of 3' 




-3 

-2 

-1 

0 

+ 1 

+2 

+ 3 

Totals 


-3 

1 

2 






3 


-2 

2 

6 

2 

1 




11 

WEIGHTS 

-1 


3 

12 

3 

2 

1 


21 

(y) 

0 


I 

4 

17 

3 

3 


28 

in S-lb 

+ 1 



1 

7 

11 

4 

2 

25 

intervals 

+ 2 




2 

S 

3 


10 


+ 3 







2 

2 

Totals 

3 

12 

19 

30 

21 

11 

4 

100 


To e\^luate 2(fxy) proceed as follows; 

(i) For each value of e.g, for x= ~3j sum the values of (j^ ) corre¬ 
sponding to it i,e. in this case (- 3) (1) + ( - 2) (2) = - 7. 

We have worked down the fimt column, and its contribution to our 
totaHs(-3)x(-7)-21. 

Work similarly for each column and put the result at the foot in an 
extra row. 

The sum of this row is the sum required. 

(ii) For each value of y, e.g. for y ^ - 2, sum the values of (fx) corre¬ 
sponding to it, i,e, here (- 3)2 + (- 2)6 + (-1)2 + (0)1. 

The contribution of this row is (- 2) x (- 20) ^40, which is put in an 
extra column to the right of the row. 

The sum of this column should also be 2 (/^) 

To determine use the values already found for 

(T^, cr^ in terms of elasi^intervais as umts^ since the F calculation has 
been done in these. 

Repeat the calculation with the same data but with the origin in the top 
left-hand square, so that all the values of y and of y are increased by 
3 and none are negative. 

It is possible as a *tour de force* to set {fx) separated from {fy)y to 
multiply by xOOOy, and to accumulate thus at the same time {fx% (fo^) 
doubled, and {fy^)* Try this out with simple figures first. 


D 


M.N.M. 
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'Irj» 2 

regarded as a set of on a weightier bar. . 

The cenm of maos, distance i from the ongtn,» g.ven by 

• S(m.-c) 

X('«)’ 

that M takes the place of/in our statistical formula, and the total 

ma^s M takes the place J about an axis perpendicular to 

“S r?.ls of Smtion about a parallel mtis through G, the 
centre of mass, we write 

/,oa.= V-{OG)*=V-^ .^ ^ 

Thus there is - e- pmaM -^ 0 “ « of 
pT Jel Ses'The pT«f of this theorem for the special case of masses m 

then /«>-/*> + M/e d/a- 

We now sum for all the values: 

2 (/k®) = ** 2/ + 2iS (A) + 2 (/»“) 

M Vf-JV- andS(f2)-0 since the a’s arc measured from the mean. 

NowX/=J''>^"“"'J ' , __„f tUp pniiation are the sums of 

fqTaS'^d S.Ji-omVorVn and the mean respectively, 

.. 
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7,8a 

Dctennine the radius of gyration of the masses in Ex. 2.9 (i) about the 
origin 

and (ii) about the centre of mass directly. 

Check by using the theorem of parallel axes. 


7,9. 

Further examples ideally suited for machine work will be found in 
statistics text-books. 

An interesting exercise is the x® testj to decide whether the character¬ 
istics of a given sample may be considered consistent with a certain 
assumed distribution of this characteristic in the ‘population* from which 
the sample is drawn. An elementary example of this test is the following: 

A die was shaken 120 times, and the distribution of scores was (1)12, 
(2) 23, (3) 19, (4) 24, (5) 26, (6) 16, 

On the assumption that the die is unbiassed, the expected value 
would be 20 for all six numbers. 

Writing for the actual and for the expected frequency for the ith 
class, X® defined by the expression 

2 ^^^^—^ summed over the classes. 

The larger the value of the less likely is the result to have arisen 
by chance for a true die,* For this case 

82 32 p 42 43 

^ 20 20 20 20 20 20 

The X® tables show that in a large number of repetitions of this 
experiment one could expect to be at least as large as this in about 
21% of them, for a true die* By common practice the assumption (true 
die, in this case) is rejected only if the calculated x® value is ‘rarer* than 
5%; and is in some doubt if it is rarer than 10%. 

[• The six ei-va!ues are not all independently variable, since their sum 
is fixed; the appropriate line of the x® table is that for ‘5 degrees of 
freedom’*] 

A further trial with the same die gives 14,19, 25,23,24,15, Does this 
make it more likely that the die is biassed ? [Solution note. 



















CHAPTER 8 


SOLUTION OF LINEAR EQUATIONS 

In this chapter we can only give the reader, by means of simple 
examples free from pitfalls, an outline of the two main lines of approach, 
viz: elimination (Exx. 8.1 to 8.5) and successive approximation 

(Exx. 8.6 to 8.9). JK A 

The numerical process of solution involves only the coefficients ana 

the lay-out consists of numbers only. We set out the skeleton data the 
coefficients and the absolute term — together with an additional sum¬ 
mation column. Thus, for the equations 

5-23x,+ 0 - 47 x 5 ,=2-63 


11-25X,+2-39x3-308 


we wnte 


(Equation E) 
(Pivotal, P) 
(Equation E') 


5-23 

11-25 

(- 0 - 0001 ) 


0-47 

2-39 

-0-6411 


2- 63 

3- 08 
1-198 


S 

8-33 

16-72 

0-577 


^3= -1-198/0-6411= -1-869 


Then fromP, Xi = (3-08 + 2-39 x l‘869)/ll-25 — 0-6708 

In line 3 we have eliminated by the following method, which is 
always adopted: 

(i) Find the numerically largest coefficient of JCj: t^ is called the 
pivot, and the equation in which it occurs is called the pivotal equation 
for this elimination. 

(ii) To eliminate beUveen the pivotal and another equation £, 
find the multiplier (M^) by which the pivotal needs to be multiplied 
before subtraction from B, 

Me is always numerically less than unity: here it is 

5-23/11-25=0-4649 to 4S 

(ill) In the next line are set the coefficients of the new equation E' 
derived from E after this process is complete, e.g, under we have 
0-47 - 2-39 x 0*4649. 

The check is that the 5 figure (obtained by the same process from the 
5-figures above it) should be a^in the sum of those in its line. 


Note. We do not here justify the number of figure retained at each stage. 
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8,1. Solve the equations 9-15x1+0*83x2 = 78-5 

0*05x1+2*46x2= -36-1 

treating the given values of the coefficients as exact: this means that the 
values of Xj, can be given finally to any desired degree of accuracy 
(say 4S). 

The S procedure checks the work as far as the equation for one vari¬ 
able; but error can arise in dividing to find this variable, or In the back- 
substitution to find the other. Substitute in the pivotal equation. 

The complete check is given by calculating to 4S the value of the 
* residual* for each equation, i,e. of 

9-15x1 + 0*83x2-78-5 
and 0-05xi + 2-46x2 + 36-1 

These should differ from zero by not more than 0*01 which represents 
the 4th significant figure in the numbers 78-5 and 36"L 

8-2, (a) Estimate the acceptance of the results for Xj and in Ex. 8.1 if 
the coefficient of the former is in fact a rounded figure, by direct 
solution with 0*045 in place of 0*05, 

A set of equations in which such a process causes large changes in the 
calculated values is called in-conditioned. An example follows, 

(b) The solutions of -+’^- = 2; 5x + ^ = 71 can be seen to be x^7, 

X v 

y = 6; and the second equation can be rewritten as =2|^. 

Consider the equations modified by rounding all the coefficients to 
2 decimal places; 

0-14X +0-17:^1^2*00 
0*17x + 0*20y = 2-37 

Solve these equations, treating the coefficients as exact. The pivotal 
multiplier for x-elimination, viz: 0T4/0-17, must be rounded (say to 5S). 
Determine the residuals. 

Notice how greatly the solutions have been affected by small changes 
in the coefficients: the equations are clearly ill-conditioned. Since in 
practice the coefficients of equations are almost always subject to error, 
the solutions of ill-conditioned equations are virtually useless. 

[Solution note 
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For o.J* ^ 

In solving a triad of equations the pivotal procedure is adopted for 
each elimination. It is vital to remember that it is always the pivotal 
equation which is multiplied, but its multiplier has one value for 
eliminating between pivotal (P) and equation E-^ and another {M^ for 
eliminating between P and £3. Since the M’s are recorded, the place to 
put Ml is against £1 and M^ against £3, although the M is in each case the 
multiplier for P. 

I'he example will make this clear. The star ^ marks the pivotal co¬ 
efficient for elimination of from the triad, and the double star ‘**'’**' 
shows the pivot for the eliimnation of from the ensuing pair of equa¬ 
tions; 


M 



X3 


s 

0*3 

-3 

2 

-0-4 

1 

-0-4 

-0*5 

5 

-4 

+ 7 

0 

8 


-10* 

1 

0-1 

0 

-8-9 

-0-4857 

(0) 

1-7 

-0-43 

1 

2-27 checks 


(0) 

-3-5** 

7-05 

0 

3-55 V 


(0) 

2-994 

1 

3-994 V 



C4S) 




=0*3340 

- 3-5) -0*6728 

= 0 * 1 x 3 - *^ 3 ) / ( -10)=0-07060 [N.B. 4S not 4D] 

Residuals (i) - 3xj 4- - 0*4x3 -1 - 0*0002 

(ii) 5 x2^-4x3+ 7^3 “ 0— -0*0002 

(iii) — lOxj^ + X2 + 0* 1X3 — 0 = 0*0002 

Solutions, to 4S 

0*07060, 0*6728, 0*3340 
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8,3. (a) Solve the triad of equations 

- 3-27x^ -r 205x2 ” 0 * 39 X 3 = P 
5 08xi - 4*13x2 + 6-95x3 = ^ 

-10*0xi + 1*00x2+ 0100x3 =r 
where p = 1, q — 0, r = 0. 

Treat the coefficients as exact, work in 5S, and give the residuals to the 
largest number of figures justified by the terms which contribute to 
them. Finally give your set of solutions to 3S. 

(b) Solve the same triad for p=0, q — l,r = 0;and 
{c) again for p— Oj q=0j r = l. 

It will be apparent that the pivots are the same for the three equations 
(so long as elimination is carried out in the same order) and the multi¬ 
pliers are the same. It would therefore be possible to carry out the 
eliminations on a single layout, with a single check-sum, as follows: 


M 

X| X, 

34 

(a) 

Cb) 

(c) 

S 

0-327 

-3-27 2-0'5 

-0-39 

1*00 

0 

0 

-0-61 

-0-508 

5-08 -4-13 

6-95 

0 

1-00 

0 

8-90 


-10-0* 1-00 

0-100 

0 

0 

1*00 

-7-90 


(0) - 

(0) — ** 

— 

— 

— 

— 

— 


(0) 

— 

— 

— 

— 

— 

Solution 

Values 


; Residuals 



of (a) 

(vii) (iv) 

(i) 

(X) 

(xi) 

(xii) 


(b) 

(viu) (v) 

(ii) 

(xiii) 

(xiv) 

(XV) 


(c) 

(ix) (vi) 

(iii) 

' (xri) 

(xvii) (xviii) 



It is convenient to insert the solutions and residuals (i?) at the foot as 
shown, in order as they are obtained. 


























NOTES FOR EXERCISES 


44 

For BA* 

It will be clear from Ex* 8*3 that if a triad with fixed ^^-coefficients is 
liable to appear for solution on a number of occasions with diJFcrcnt 
absolute terms, then the grid of solutions obtained by your work will 
provide a basis for a calculation of any case. If we write the solutions of 
(a) as Xya on, the solution of/»= 7, ^ = 5, r = 0 is 

Xi= = 7x^ + 5x^ ; =7X3^ + Sx^j, 

because - 3 -27xj^a + 2*05x2^ - 0*39x^o = 1 

- 3 -27x13 + 2*05x33 - 0*39x33 = ^ 

- 3*27(7xi, + SXi 3 ) + 2*05(7x2, + ^^ 2 .) -0-39(7x3, + 5 x 33 ) -7 

and so on for the other equations* 

A grid of nine figures like that forming the L.H.S, coefficients of 8.3 
is called a 3 x 3 matrix. The R.H.S* coefficients also form a 3 x 3 matrix 
w^hich is here a ‘unit matrix'. 

The solutions also form a 3 x 3 matrix. This aspect of the problem is 
considered again in Ex. 8.5. 

After 8,4. Non-regularity for a triad. 

Failure to obtain a unique solution can clearly occur if the final 
elimination gives 0, Xg d* 

If d=0 (my value of Xg will satisfy, and if no value. That there is 
failure depends on the left-hand coefficients; but the type of failure on 
the right-hand also. 

For a triad the situation can be considered in terms of coordinate 
geometry in three dimensions, in which an equation of the first degree 
represents a plane. 

Immediate failure occurs if the two first equations give planes parallel 
(no solution) or coincident (an unlimited number of solutions). 

Failure at the end of the calculation occurs if the line of intersection 
of the first two is parallel to the third (no solution) or lies in it (an 
unlimited number of solutions). 

For further reading. 

(See Bibliography for full titles.) 

Matrices; introduction and manipulation: 

Matthews, Matrices (I) and (11) 

Matrices in the general mathematical setting; 

Moakes, Core of Mathematics 
Linear equations: 

Cohn, Linear Equations 
Numerical methods: 

Wooldridge, Computing* 
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8.4. Deduce from your triple solution of Ex. 8.3 the solutions of the 
same equations for the case p = 5*20j q = -1*74, r 3*87, Use these 
values as multipliers for the solutions obtained with unit right-hand 
sides, as shown in the note. Calculate the residuals. 


8.5- Solve the two pairs of equations, using a single layout 

523x1 + 047x2-1*00 -0 

11*25x1+2*35x2=0 -1*00 


Deduce the solution of the pair already solved on p, 40, by the method 
of multipliers as above. 

This method is of great practical importance, and we shall set it out 
(for two equations only) using the algebraic notation with suffixes. 

The square array (*matrix’) of left-hand-side coefficients, which 
determines the solution procedure, will be written 

/ <?ii 

\^21 <^22/ 


The first suffix in each case refers to the equation, upper or lower; and 
the second suffix to the variable, Xj or Xg. 

The solutions of the equations, solved first with and then with 

as right-hand side, will also be written in square matrix form: 

/«ii “ia\ 

\«21 ^22/ 

Here the first suffix refers to the variable, and the second to the solution 
— first of then of ^ as right-hand side. 

By substitution into the equations we obtain: 


^ 12^21 ^ 11^^12 + ^ 12“22 ^ 

a2i«ii + flggagi = 0 J + ^32^22 ^ ^ / 

These are the conditions that the product of the (2-matrix into the 

a-matrix is the unit-matrix ; i.e. that the a- and a-matrices are 

inverse. The process of obtaining the a's is called inverting the matrix of 
coefficients* We have seen in Ex, 8.4 that the r^ult enables us to derive 
the solution of any set of equations with these coefficients, whatever the 
right-hand sides may be. 

(See references on left-hand page, for further reading.) 










46 

For 8.6. 


NOTES FOR EXERCISES 


Solving Linear equations by successive approximations. 

We have seen that the acceptability of solutions can only be shown by 
slightly varying them and watching the effect on the residuals. The 
following method follows this approach, starting with arbitrary values of 
the variables. 

Consider the equations 2x + 5y-8=R.; 

30x +y i 28 

where we require y such that = 0. 

We start with a schema, in which the coefficients figure in a new role 
(and rearranged^ viz: to determine the changes made by 

changes ^x — l,Ay — l in the values of the variables. As a check we 
calculate at every stage also the value of Our schema is 

AR^ AR.^ AR 

(a) Ax = l 2 30 32 

(b) Ay = l 5 1 6 

Note that line (a) contains the first column of coefficients, and line (b) 
the second colunm, from the given equations. 

We now start with, say, .v — 0 andy = 0 and trace the solution through: 
at each stage we reduce one of the residuals to zero or conveniently near 


zero. 





R 

Initially 

* = 0 

Ax 

y=0 

Ay 

-8* 

28 

20 

By (b) 


1-6 

0 

29-6* 

29-6 V 

(a) 

-0-99 


-1-98* 

-0-1 

-2-08y 

(b) 


0-4 

0-02 

0-3* 

0-32^ 

(a) 

Solution 

-0-01 
—1-00 

2-00 

0 

0 

0 7 


* The starred figure in each case controls the next move. 

The operative values of Ax^ Ay have been chosen as follows: 

(i) To add S, to make zero, Ay = S/5. 

(ii) To reduce- 29*6 to zero exactly we should take A x={- 29- 6)/ 30, 
but we prefer at this stage the 2D value Ax— - 0'99. 

(iii) Again, as in (ii), we shall not allow a complicated value of a 
change to be used. Take Ay = 0-4 exactly^ to bring /?j nearly to 
?:ero. 
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8.6, Repeat the example given, making as the first step a change in x to 
reduce R^^ This is in accordance with the working rule, not followed in 
the example^ of reducing first the largest residual. Notice that it is not 
necessary to reduce to Progress is made, without introducing 
heavy figures at an early stage, if we can reduce the magnitude of by a 
factor of about 10: thus here it would be quite enough to reduce R^ from 
28to -2. 

The strict method of solution, in which the R'^ are made to vanish 
alternately, can be seen by plotting graphs of the loci to be ^spiralling* 
round the point of intersection if the gradients are of opposite sign, 
switching from one line to the other and back: 


x = OA^ 

y = 0 

is a pair of trial values. 

OA^ AB is the next 
pair. The coords of C, 
the next; then of /), 
and so on. 

When made to vanish we obtain a pair of coordinates wliich 

satisfy line 1 only: one coordinate is carried over to obtain a point on 
line 2, and so on. 

If, however the residuals are only made approximately zero, one is 
working with a band rather than a Hne. In the end stages the residual 
must be made zero. 

It will have been seen that the speed ofthis method of solution depends 
a great deal on judgment in choosing which residual to reduce fii^t and 
by how much. There is nothing to prevent the solver from building up 
modifications in any way — including a change of two or more values at 
once — and freed from its mechanical character this is known as the 
method of Relaxation. See Note 8.8. 

For computer work, on the other hand, the purely mechanical routine 
of the spiral is appropriate; but the order of operations must be chosen 
in such a way as to ensure that the solution converges, (It may do so by 
steps, not by a spiral.) 

The same method is available when one or more of the equations is 
non-linear: in fact there may be no alternative in such cases. 

For further reading see Wooldridge pp. 115 seq. 
























4g NOTES FOR EXERCISES 

For 8.7. If the absolute term had been 28-1 in place of 28, the columns 
and R would have been different, and after the fourth Ime (fJ*=04* 
/i =042) we should multiply temporarily by 1000 and proceed as follows 


By (a) 

(b) 


400 

20 

400* 


420 

-13 

-6* 

10 


4 7 

1-2 

0 

11*2* 


11-2 7 

-0-37 

-0-74* 

0*1 


-0-647 

0-15 

0-01 

0-25 


0-267 

-0-01 

-0-01 

-0-05 


-0-067 

e relaxation scheme 


dRi 

JRs dRj 

dR 

i 

-S 

2 

■3 

-6 

Ay^i 

6 

7 

8 

21 

Az = \ 

-1 

4 

2 

5 


(a) 

For 8.8. 


W 

before going on we form combinations which will relax two or if 
possible only one residual at a time, e.g. 

(d) dx = l My = l 
(c) Jy=-; 

if) ^*=2 

Thenifthedat: 

3c=0 y=0 

By/; Ax= - 2 
By e: Ay^ - 


Az=l 

0 

13 

7 

207 

II 

-10 

9 

0 

-V 

Az= -1 

-9 

0 

-8 

-177 

^ve, for: 

Ri 


Ri 

R 


1 

-9 

-8 

-16 

Az=+l 

10 

-9 

0 

17 

.,Az=^ 

0 

0 

0 

07 

is 





-2,y = 

-1,^ 

= 5. 




For 8.9. o . c 1 

If however the residuals Rj etc. above were 1 , — 8, — 8 the final step 

would give 0 1 u- 

We have too residuals zero, but no combination which reduces this. 
We obtain it thus; combine d and / to get another set ending in zero 
(it is 8d + ?/). 

(g) Jx=22,dy = 8, = l -63 104 0 41 ^ 

This minus (e) X 6*3 is 0 47-3 0 47*37 

i,e, Ax =22, ^ — 24'2, 

and this final modification x { - 1/47-3) gives the required solution. 
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8p 7, Verify that the mechanical process on the left hand page gives 
5c=-1 00338 ^^2 00135 and show by direct substitution that the 
residue given are correct. Repeat Ex. 8.6 with 28^1 for 28 and continue, 
to get a solution with residuals less than 10^. 

8 .8 p 

(i) State the triad of equations for which the given scheme is 
appropriate, and check directly that the solution satisfies them. 

[Solution note 

(H) Solve the equations with the same jc-coefficients but for which 
the residuals (for zero jc, are -10, -4, -7* 

(iii) Solve the equations 

K+y+2 =09 
10i£-y +132=1*0 
7x + 2y-32=3*8 


8*9, 

(i) Solve, to the accuracy required to obtain residuals numerically 
less than 10”® 


-Sxj+dXfl-Xa =10 

2X|+7x2+4x^^ 5 
- 3xj + 8 x 3 + 2 X 3 = 8 

(ii) Solve, with residuals numerically less than 10“®: 

X 1 +X 3 +X 3 = 1*0 

lOxj -X2 + 13x3 = 1*0 

7xi + 2 X 3 -3x3=4*1 
















CHAPTER 9 


NUMERICAL INTEGRATION 

This process is equivalent to finding the area ‘under’ a curve y =/(■*) 
between given ordinates, by calculation of a series of ordinates at inter¬ 
vals throughout the range; e.g. to find J f[x) dx by dividing the range 0 

to a of * into a number of parts, usually 2» of width k so that a=2nh. 

Elementary methods are (i) the trapezium rule, by which each area of 
width 2h is treated as a trapezium. The first will have area 

K/(0)+/(2/.)}2A; 

and so on. (ii) Alternatively we could take the mid-ordinate as the 
effective height of each 2/j-strip; and give the approximate area as 
/(/j) X 2A for the first and so on through the ordinates at odd multiples of 
A. 

A little examination of curves will show that when one result is too 
small the other is too large and vice versa, and the method most used in 
practice combines the two. 


Simpson’s Rule states that a good approximation to the area of the 
first 2A-strip is {^/^K/C^) + 

This is equivalent to writing 

- J{/(0)+ wo +m)) . (s) 


Similarly 

j2mdx=^^{f{2h) +4/{3A) +/(4A)} 


and so on. 

Finally by addition 



dx 


h r/(0) + 4{/(A) +/(3A) ++/(2rA ^ A) +... ^f{2nk - A)} 1 

”3 L + 2{/(2A) + /(4A) ++/(2rA) +... -^f{2nh -2h)} +/(2nA)J 

where 2nh = a; and the result is more accurate the larger is the value of n 
and consequently the smaller the value of A. 
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9.1. 

Determine by Simpson^s rule an approximate value of 


^ 

(ji+x 


taking h = O'l, Employ the layout below: 


X 

1-fX 

f(x) values 

multipliers 

U 

1 

1-000,00 

1 

0-1 

M 

0-909,09 

4 

0-2 

1-2 

0-833,33 

2 

0-3 

1-3 

0-769,23 

4 



etc., 

to: 

1-0 

2-0 

0-500,00 

1 


(Space for differ¬ 
ences if used) 


Determine the cumulative sum checking that the cycle register 
shows the total 30. Why ? 

Then the value required = A2/3 

Compare the result, rounded to 4S, \nth the true value of log* 2 

9.2. Repeat the above calculation with h = 0*05 

9.3. With A =0‘1 find a value for 

ri dx 

I i - ^to 48^ working to 6D 

and compare the result with the knowTi value of 7r/4 to 4S 
Note: the routine, of carrying tv\^o places more than are finally required, 
does not necessarily secure an accurate result when an approximate 
formula is used. The adequacy of the method, with any correcting 
devices included, has also to be considered. 

dx 

9.4. Find a value of -—~ , working to 5D 

J Q 1 + 

(i) with h = 0*2 and (ii) with h=0*1. 

It will be advisable to form the third column with all values off (x) for 
intervals 0*/, and to extract the required values for the two calculations 
separately. ’Fhe f(x) values will be needed again (Ex. 9,6). 

The comparison of the two results will give an indication of how many 
figures in the result can be relied on. 

Compare the result with (1/3) log^ 2 + 27r^3/9. 
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For 9.5* 

It will be seen that if the integral is one which can only be evaluated 
numerically —- and it is for these that the method is most useful — one 
might have to face a series of re determinations in order to get the re¬ 
quired accuracy. 

This is the method used on an electronic computer, but for desk work 
we must know for any chosen value of h the order of error to be expected, 
and if possible a correction which we can apply. 

We start with the fact that the basic equation (S) of p, 50 is exact if 
/(x) is a polynomial inxoi degree 3 or l<^s. In verifying this it is easier to 
t^e the origin at the centre of the range of x, so that we prove: 

/(*) =*{/(- A) + 4/(0) +/(A)) 

The reader should write f{x) = establish 

that there is a discrepancy proportional to This is zero if p is zero, 
and can otherwise be reduced to acceptable size by diminishing the 
interval h. 

Looking at the Simpson method for a range a = 2nh^ wt see that if f(x) 
can be represented well enough by polynomials of degree 4 — they can 
be different polynomials for each range 2h —■ then the correction to be 
applied is 

where p is the mean of the values of jp. 

We have seen that the 4th differences for such a polynomial 
function will be {4l)ph^. We may therefore rewrite the correction as 

C = — 

where ^ is the mean fourth difference, 
since a = 2nh. 

A value for ^ can be assigned sufficiently well from the (2n - 4) 
values of the 4th differences which will arise from the calculated values 
of/; and which will have been determined already in order (i) to check 
the accuracy of the/ values; and (ii) to establish the adequacy of this 
method, for the chosen value of A, 


CORRECTING RESULTS OF SIMPSON^S RULE S3 

9*5* 

Form a difference table for Ex, 9.1 and compute the correction: 
compare the corrected value with the true value. 

Estimate, without recourse to a new table, the correction to be applied 
in Ex, 9.2. The 4th differences, being of the order of will be approxL 
mately of the previous values; hut there will be twice as many, 

9*6* Tabulate the function of Ex. 9.3 and apply the correction calculated 
from 4th differences. 

Notice that for small x, 1/(1-hx^) 1so that the earlier ziy 

should be 24/r^; but this expansion is only appropriate near 
Near x = l write = 1 - ^ and we have 

= iO-h(z-x^/2) + (z~z^/2)^ + ( y + ( 

Examination of the term in shows differences to be 6A^ in this part 
of the table, 

9.7. Use the tabulated values of Ex. 9.4 to form a table of 

f(x) = l/VCl+x3) 

for the same range. 

Use an iterative method for this square root starting with the value 
from 4-hgure tables; or as in note below. 

Form 4th differences, correct your Simpson result and indicate to 
how many figures you consider your result reliable. 

Your check can be eMer two corrected Simpson results for different 
/t; or a different integration formula. 

AW.' Working directly from values of A"=(l -hx^) you could take one 
move only to find Write F(x) = (1/x^) 

Verify that Newton *s rule gives 

»«-n=y(3 ~iVV) 


E 


M.N.M, 




















CHAPTER 10 


NUMERICAL DIFFERENTIATION 

Consideration of the graph of a function shows that integration will 
give a result to rather % accuracy than the w^orst of the data: e.g, 

if a function is represented to 1% or less over the w^hole range, then the 
Simpson formula will generally give a result well within 1%* 

On the other hand, differentation gives a result with smaller accuracy. 
This will be familiar to a physicist who has drawn a cooling curve for a 
body by plotting temperature against time, and has tried to derive from 
it a series of values of the rate of fall of temperature. He finds that con¬ 
cordant results can only be found by ‘smoothing^ the data, A chord 
drawn to join successive points will indicate by its gradient a good 
approximation to the rate of fall at the middle of the time-interval, 
provided the data are knowm to a high enough accuracy to give a smooth 
curve. 

The criterion of smoothness for a series of numbers is the regularity 
of the differences. However, nth differences can always be irregular to 
the extent of (1/2) x 2", since the final figures of the data are uncertain to 
the extent ±(I/2) ; differentiation therefore requires say 3 extra figures 
in the data in place of the usual 2, to get a result to a required number of 
figures. 

We use as a first approximation 

or nx+hjl) + -/W) = (1/A)J/W 

If a difference is regarded as belonging to the value of ike Junction 
opposite which it is wriitm the symbol 3 is used, and we write 

f\x) =- {\ik)hf{x) ...,{10.A) 

Notice that if / is tabulated for ^= 0 , 0 - 2 , 04 etc, up to 1-0 the 
differences He opposite the missing values 0-1 etc. and the formula gives 
(approximately) /'(0ri),/'{0'3) to//0^9). 

Consideration of polynomials gives the improved formula 

fix) ^ (l/A){S/~ (1/24) 8^+ (3/640) S^}.(lO.B) 

in which omitting the last term is less serious usually than errors due to 
the data, the values of/. 
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10.1. Referring to/(x:) — ^ 


tabulated and differenced in Ex. 9.5 for 


0(0'2) I'O, obtain values oif{x) for 0-3 (0'2) 0-7 

(i) from a graph, using 3S for/(%); 

(ii) using formula (10,A) or the simpler forms given above it; 

(iu) using (10,B)j remembering that h^f means the value on the 

same line as the value of being considered; 
and (iv) by calculating/^ (x)= - 1/(1 +a')^. 


10,2. If f(x) = + ajX + -s- agX^ + a,jX^. 

Show^ that f [k) -/( - ^) “ S/(0) = 

(We are taking 2k = /i, the tabular interval) 

(iz) Hence show’^ that the error in (10,A) is or 
(b) Find S/(^) =/(2A)-/(O) 1 and hence evaluate 

and fi/(-^)=/(0)-/{-^) / Sy(0) = a/(A)-S/(-A) 

Show^ that /"(0)=^(l/A^)Sy(0) with an error of order 
(Notice that k—hl2 is only used as a matter of notation: no 
intermediate tabular values in the A-intervals are used), 

(^r) Using 4-figure x^lucs of sin x over 0®(10^) 90®, form a difference 
table. Obtain values of/" (,v) for x in radians by w riting k =ir/18 in 
the equation 

f"{x) (l/li^)S=T(x)„.,. .. „(10.C) 

and compare with values of “/(a:). 


10.3. Form a table of x to 6D, x = 1 (1)10. 

Difference and find U(x) for 2'5{1) S'5 by (10,B) and compare with 
values of (4343)/x. 
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NOTES FOR EXERCISES 


For 10*4, This example was important because if/ could be represented 
by a formula it would be possible to predict with some confidence the 
behaviour of the substance beyond the range of the experiments. The 
design of power plant depends on this function for steam. 

For 10,5, 

Formula (10.C) can be rewritten in the form 

hT (*) ^ ={fix + h) -/(*)} - {/(*) -fix - h)} 

=fix+h)-2fix)+fix-h) . (10.D) 

Thus if we are given a rule by which/''( aj) can be evaluated for any x 
(lO.D) giv^ a method by which the values of/(jS”A) and f(x) can be 
used to find a value of f(x + k)y with an error of order This is one of 
several methods for numerical solution of differential equations of the 
second order, (There are devices by which the correction terms of order 
h* are fed in; either on the way, or en masse after tabulation of the 
approximate solution), e.g. if 

r(x)= -fix) 

then 10(D) becomes 

fix+k)^i2-h^)fix)-fix-h) 

Layout as follows: 

X x+h f(x-h) f{x) /(jc+Zf) approx, 

(1) given; and (2) givens (3) calculated 

(2) given (3) from above,—^ (4) 

and so on. 

The full set of values, after the given pair, reads down the final 
column, the values of the variable being those in the second column, 

(The arrows represent work which in some cases might require a 
further column in the layout.) 
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10.4. *Live’ examples of numerical differentiation are rather rare. The 
following is one from physical chemistry: 

The energy change involved in the dissociation of a molecule of a 
certain compound is a function /{T} of the temperature (^*0) at which 
dissociation occurs. This function is required, and can only be found by 
determining experimentally another function ^(T) which is known to be 
connected with it by the equation 

dT {T + inf 

You are given the smoothed experimental data as follows; 

T 5 10 15 20 25 30 35 40 45 50 

q 1*700 1'962 2-175 2-347 2481 2-581 2-652 2-700 2-733 2-767 

Taking A = 10“^, form a table of/for 3^=12- 5(5) 42-5 and by differen¬ 
cing this table state whether/ could be representeeJ by a polynomial over 
this range of values, [S olution n ote, 

10.5. Solve approximately the equation f*{x) = -f{x) to find values of 
f{x) at intervals of 0-1 in x given the initial vaiues/{0) = 0 and 

/(_0-l)=-0-099833 

using the relation, applicable to this equation only 
fix+k)^i2-hYix)-fix-k) 

Errors arc accumulating, and it will be possible to see at what point 
the value of/(:v) differs from sin x in the 4th place. 

10.6- Using the expansion 

cos (h) ^ 1 hS/2 + hV24 - h®/720 

show that, to order 

cos {x + A) = cos x{2 - A“ + A^/12} - cos (x - A) 

and hence repeat the process of £x. 10,5 with the correction ‘built in*, in 
the revised multiplier. 

Start with/(O) = 1 and/(A)^/(0-l) evaluated by the expansion. Work 
with all the figures your machine SHOW'S, up to cos (1-5) =0-070737 to 6D, 
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PROGRAMMING AND THE COMPUTER 

1 * In these exercises the art of programmmg a computation has been 
developed. It is not sufficient to have a method which is sound; broken 
down into ‘units of process* which the machine can handle^ it should also 

(i) proceed smoothly and economically; 

(ii) lend itself to a clear layout which makes errors less likely; 

and (iii) be capable, if long, of intermediate checking, to avoid wasted 
labour, (Study the example in 6,2,) 

In non-iterative processes, e,g, solution of simultaneous equations by 
Gauss’ elimination method, a checking method is built into the process 
itself. 

Iterative methods are self-checking; and the worst that can happen 
after, say, a copying or re-entering error is that more iterations are 
required before the result is achieved. 

It ^vill be noticed that the programmes which lend themselves to 
flow-diagram representation are all iterative* The instructions which 
form part of a loop are in fact repeated, with modifications which are 
provided for, until a certain routine test is satisfied and ‘directs’ the 
process out of the loop. Look again at the diagram on p, 16 from this 
point of view% 

Electric desk-machines 

An electrically-driven machine which provides for automatic division 
is in fact designed according to the Figure on p, 16, No bell actually 
rings, but a close scrutiny will show that every subtraction actually 
continues until zero is passed, and the last backward turn is then 
retraced* 

Some electrical machines will only multiply digit by digit: pressing 
the ‘4’ key, for example, simply causes the crank to turn 4 times. On 
fully automatic machines an extra register holds the multiplier digits, 
preselected: the ‘ x * key then initiates multiplication as an iterative 
process, continued digit by digit so long as there arc non-zero digits 
shown in this register. 

These machines are not in any true sense computers; they simply 
carry out one, or at most two, planned patterns of calculation. They do 
however suggest to the mind the desirability of machines which could 
carry out not only more flow-patterns but preferably any required type of 
pattern, according to a programme of instructions fed into it in advancet 
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Around 1840, Charles Babbage worked out all the main principles of 
such a machine, to be fed with instructions on punched cards, and the 
partly constructed machine (purely mechanical) is to be seen at the 
Science Museum, London, It is possible today for a young student who 
has some experience of calculation with a desk-machine to list for 
himself some of the basic features he would wish such a machine to 
have, if the intervention of the operator is not to be required during the 
programme. 

Requiremefitsfor a Computer 

Babbage called the part of his machine where numbers are currently 
being processed, the MILL, It was required (like the automatic desk- 
machine of today) to add, subtract, multiply or divide pre-set numbers. 
We also require it to compare two numbers: the test for leaving a loop is 
commonly dependent on such a comparison (e.g* on p, 18, ‘Is T 
negligible?’ in precise terms might be Ts | T [< 10'*?’), 

In place of our pencil-and-paper we need STORAGE facilities; for 
data, for intermediate and final results, and (not least) for the actual 
programme itself and any ancillaries associated with it* 

Above all, to obviate the interference of the operator at each new 
phase of the programme we need a CONTROL UNIT* Two factors 
have contributed jointly to the success of this, the characteristic feature 
of the computer, viz: (1) algebraic logic, initiated by George Boole in 
1840, and (2) the invention of electronic pulse-operated devices* 

The relation betu^een the parts of the computer and its input and 
output units can be sketched thus: 


INPUT OUTPUT 



One can use a computer successfully with only slight knowledge of 
the inside processes, by storing in the computer a suitable ‘translator 
programme’. This allows the programmer to feed into the machine a 
set of coded instructions not far removed in essence from a flow diagram. 
While being read in, these arc checked for translatability into machine 
process; and when the whole programme is stored it aw’aits the final 
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button to set it into operation. (Compare how, in an electric machine, 
the feeding in of a multiplier does not itself cause the machine to 
multiply.) 

The student who wishes to get an introduction to computer program¬ 
ming will be well advised to go to a short course at a college which has 
its own computer. The choice will be either a course related to a particular 
range of computers (e.g. Elliott Autocode for the Elliott rangCj or MAC 
for LC,T, computer range) or a course in one of the generalised 
computer languages (Algol, or Fortran). Which type of course is prefer¬ 
able depends on the object in view. 


APPENDIX II 


INSTRUCTORS’ NOTES 

STANDING ORDERS: for students using the machine. 

(1) See that the machine is safely placed, i.e. on a firm fiat table or 
desk. 

(2) See that you are in no danger of being jogged while at work: no 
other person should be overlooking your work unless he is 

either teaching you 

or being taught by you. 

(3) Always complete one movement before starting another. 

(4) NEVER FORCE ANY LEVER: in case of jamming — put 

machine av%^y and report the 
difficulty, 

(In many instances what appears to be a jam is only a stoppage due 
to etVAfr failure to complete a turn or failure to effect a full 
clearance. The remedy in each case is obvious,) 

(5) After use: clear all registers 

return carriage and "1 , 

locking levers ) to normal, 

and put the machine away unless directed othen^tisc. 

Further notes for class-imtruction 

If you use a machine for class-teaching with junior students it is 
advisable to rehearse in detail the calculations proposed. 

The students will naturally be anxious to carry out manual operations 
themselves, and it is tempting to devise a cooperative calculation such as 
adding the ages of the students in months, letting each one add his own. 
In fact such a project is very liable to go wrong, unless a previously 
briefed student is standing over the machine checking that every result is 
entered correctly. Individual projects (e.g, for each to turn tus age in 
years and months into months) are better because (i) one operator's error 
doesn’t affect a final result and (ii) the worker will check the result on 
paper and will be able to compare his speeds with and without the 
machine. It is still a good idea to have a student supervisor to minimise 
error. With a lock it may be practicable to set the same figures in the 
locked part, so that a running total also is obtained, 
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MACHINES 

A number of companies distributing in the U.K. have hand-machines 
with full tens-transmission (essefitml) and back-transfer (demabie) at 
educational prices in 1965 from £35 up to £58. These all have lever 
setting, keyboards being more expensive. Their capacity is usually 
10 (set) by 8 (cycle) by 13 (accu,). Choice will depend not only on cost, 
of purchase or hire or hire-purchase, but also on the service facilities 
available in your area. 

Some suitable models are listed below; the order is alphabetical (of 
distributors* names): 

Multo: Addo Ltd,, 5 Worship St,, London E.C. 2* 

Teacher 11: ADM Ltd,, Chandos House, Buckingham Gate, L 

Odhner: Block and Anderson, Cambridge Grove, W. 6. 

Minerva Master: Broughton*s, 6 Priory Rd*, Bristol 8, 

Mentor: Muldivo Ltd., 4 Dorset Bldgs,, Fleet St,, E,C, 4, 

Brunsviga 13; Olympia Ltd,, 35 Red Lion Square, W,C. L 

Keyboard {including ekctrical) machines are sold by some of the above 
and also by: 

Monroe Calc. Mach, Co,, Litton Ho,, Aldersgate E,C, 1 

and Olivetti Ltd,, 29 Berkeley Sq., W, L 

A pocket machine 8x6x11 (no transfer) is the Curia I: London Office 
Machines, 5 Lower Belgrave St., S.W. 1, 

(All addresses are in London unless otherwise stated.) 
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SOLUTION NOTES 


1.6 Turn once in first position and once in second. 

1.8 First difference = + 3w +1) + b{ln + !)+£; 

second = a{bn + 3) + 26; third = ; fourth = 0, 

3.5 (i) km, total x O'OOl =0'10, to 2D. 

(ii) max, error 0-25 km. =0*16 miles, 

3.8 (i) 0*895 ±0*009 

(ii) 5*25 ± 0*01 since x^y uncertain to 0*005 

(iii) xAy -k-y Ax = 0*005(x + 7 ) =0*026 to 2S 
.% say xy ^ 3*300 ± 0*026 

(iv) 2xy Ax + x^Ay^ 0"005ii-(jc + 2y) = 0* 135 
giving = 14*916 ± 0* 135 

Extreme values directly: 15*050,14*779 

4.1 S„ = l+r±r^ + ...{n terms) and but a 

more useful relation is = 1 -f rS^ which gives a simple 
programme of calculation starting with Sq = 0 or = 1 . 

4.3 30«-29« = (30-29)(30 + 29)-30 + 29 

and so on: thus we are summing the integers 30,29., .21, 


4.6 (a) 


^m= 


(« + l)W 
rl 


r! ’ 






(r-l)l 


4.8 (i) Aa = 0*005 for every a, and thus 

errorii^0*005(x^+x^-h... + 1) =0*005 x8*21 

= 0*041 to2S 

(ii) error in ^ 4x® dx, and so on. 

Calculate (4flQX^ + 3a|X^ + la^x + a^) Ax 
or recalculate/(x) for x = 1*26, 

4.10 The settled value of 4th difference seems to be 24, and the errors 
presumed in the other values ( + 9, -36, +54, -36, +9) 
indicate an error + 9 at the ‘centre of disturbance*. The reading 
6*10 should be 6*01: transposition of figures in copying was a 
likely cause. 
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5.1 (b) The rounding error for 4S is anything from 0 to 5 (or - ) 
in the fifth figure. Taking x as exact, a: x (l/jr rounded) could be 
in error by as much as 1*249 x 5 x 10"^ for the first result and 
0*08334x5 X 10"^ for the second: thus rounding is potentially 
seven times as great a source of error in the second case as in the 
first. 


5.2 Max. error in numerator xy due to rounding 

= 39*2x0*005 +2*47 x0*05 
=0*320 

in yyjz ^ 0*320/43 = 0*007 

Cf. greatest and least values by direct calculation. 

5.3 (ii) A general result for which is often quoted in 

scientific work is showm as follows: 

We have log^ z = log^ log, x-\-n log,^. 

Then if x alone suffers a small change Axj we can use the result 

\ (iz_m 
z dx X 


to obtain 


Az m Ax 

Z ~ X 


[Fractional or % error in z^ due to x alone, is | r/t | times the 
fractional or % error in x,] 

A similar result is true for y alone; and if we neglect their 
interaction we obtain for the greatest numerical error in z the 
equation 


max 



Ax I 

m “ ^ 

x j 


I Ay 

n — 

\ y 


We have seen this in the special cases sr=xy, z = xjy. 


6.1 It is worth while to consider whether build-up division gives an 
easier programme, [Set a*; accumulator clear and fully to right; 
crank and shift so as to get N in accumulator correct to as many 
figures as possible,] In cither case the flow diagram for is 
as follou's: 
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This flow diagram can be given more precision by stating (pro¬ 
gressively) the number of significant figures to retain at each stage. 


7.9 Yes. Although the second trial by itself gives a smaller value of x\ 
the two regarded as a single longer trial gives a much larger value; 
the persistently low frequencies 1 and 6 contribute largely to this, 
and the result accords ^^^th a common-sense view. 


8.2 (b) Any lingering doubts about the sources of error can be 
dispelled by an ’elementary arithmetic* solution with exact 
multipliers 14 and 17, giving x = 3t, . 

Graphically we see that the lines are nearly parallel, so that a 
slight change has a large effect on the point of intersection. 

Algebraically, the determinant 


with the actual values of the a*s. 


■u 

a^i 


Hz 


is small compared 


8.8 (i) 


10,4 


dT 


-Sx^6y-z^l=0t etc. 

requires knowledge of and is therefore only obtainable 


over the range stated. Even so, its acceptance is scarcely 3S 
(below 32'5) or 2S (above); and the same is true of/(T). We do, 
however, obtain a steady value of for this range and to this 
accuracy: thus / can be expressed well enough by a cubic 
polynomial, either in terms of T or (as the chemist prefers) in 
terms of absolute temperature T'+273, wTitten below as 

The above w'ork indicated the need for a greater range of 
experimental values and also an improved numerical method. 
All this has now been achieved, as below. 
q was fitted to a formula of this type: 




(This corresponds to a qtiadralie polynomial for/) 

The method was (i) to compute and log^,^) 

over the range; (ii) to find, graphically, the value of such that 
g — aik is sensibly constant; (iii) knowing that $q-ai0 log^ B is a 
quadratic function, vix: fit the values of these 

coefficients to three temperatures spread over the range of 
values; (iv) check the fit of the q formula over the whole range, 
and hence confirm the formula for/= - \- aiB + a^B^. 

This formula agreed better with results obtained in other ways 
than did a cubic of slightly closer fit. This indicates that one 
would have be^n justified in smoothing the experimental data 
more than ’was actually done. 
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(Reference is to example-numbers unless a page*numher is stated,) 


Acceptance, 3.8 
Accumulator, p. ix 
Accumulation of quotients, 5.4 
Adequacy, 9.3 

Binomial expansion, 6.S 
Brunsviga, pp, ix; 3.5 

Centre, of mass or gravity, 2,9 
Checking, p, 22 
by differences, 1.7; 
en route, p. 58 
Classification of data, 7.1 
Clearance, p. x; p. 58 
Coinage, obsolescent British, 3.2 
Compound interest, 3,6 
Complement, t.4 
Continuation, 3.1 
Control unit, 6.2, p. SB 
Correlation coefficient, 7.7 
Cosine, see Sine 
Cycle register, p, ix 

D (abbreiiiation), 3.4 
Decimal point marker, 1.2 
Differences, 1.2; 4.4 
correction by; 

for integration, 9.5 
for polynomial, 4.10 
Differentiation, numerical; 10.1 
Division, build-up; S.4; tear-down, 5.1 

e, calculation of, 5.7 
Elimination, (Gauss*s method), 8.1 
Error, greatest possible, 3.8; 4,8 
detection of tabular, 4.10 

Flow diagram, pp. 11,16,18, 64 
Frequency (table), 2.5; 6.2 

Gauss, see Elimination 
Gauss-Seidel, see Iterative method for 
linear equations 
Gregory-Newton formula, 4.4 
Gyration, radius of, 7.8 


Ill-conditioned equations, 8.2 
Inertia, moment of, 7*8 
Input, p, ix; p, 58 
Interpolation, 4,4; 6.7 
Iteration, for linear equations, 8.6 
Newton’s formula for, 6.2 
for reciprocal, 6.3 

Jamming, p. x, p. 60 

Lock, p. X, 2.S 

Logarithm, calculation of, 5.9 

Matrix, 8.5 
Mean value, 7.1 
Mean, *falsc’, 7.3 
Metric equivalents, 3.3; 3*4 
Mid-ordinate rule, 9.1 
Multiplication, 1.6; 2.1 
Multiplier register, p. ix 

Operating handle, p. ix 

evaluation of, 5* 11 
Percentage, 3.1; 5.4 
Pivot; pivotal equation, 3.1 
Polynomial, evaluation of, 4.8 
building up values of, 1.7 
Present value, 4.7 
Product register, p, ix 
Product, triple, 2.2 
Product-moment, 7.7 
Programme, 6,2; p. 58 

Recurrence formula, 4*1 
Registers, p* ix 
Regularity, 8.1; 8.5 
Relaxation, 8.6 
Residuals, 8*1 
Rounding, 3*5 

B (abbreviation), 3.4 
Sequence, 3*1 
Setting register, p. ix 
Shift, p* ix 
Simpson’s rule, 9.1 
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Simultaneous equations, 8,1 
Sine, evaluation by series, S40 
by step-by-step process, 10,6 
Standard deviation, 7.3 
Storage, 6.2 :p. 58 
SufHx notations, 8.4 
Summation of some series, 4,1 
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Transfer, p, x 
Trapezium rule, 9.1 

Units, conversion of, 3.3 

Variance, 9.3 
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